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ABSTRACT 

The study of the planet-debris disk connection can shed light on the formation and evolution of 
planetary systems and may help “predict” the presence of planets around stars with certain disk 
characteristic s. In preliminary analyses of sub samples of the Herschel DEBRIS and DUNES surveys, 
Wyatt et al. (20121 and Marshall et al. (2014) identified a tentative correlation between debris and the 
presence of low-mass planets. Here we use the cleanest possible sample out of these Herschel surveys 
to assess the presence of such a correlation, discarding stars without known ages, with ages < 1 Gyr 
and with binary companions <100 AU to rule out possible correlations due to effects other than planet 
presence. In our resulting subsample of 204 FGK stars, we do not find evidence that debris disks are 
more common or more dusty around stars harboring high-mass or low-mass planets compared to a 
control sample without identified planets. There is no evidence either that the characteristic dust 
temperature of the debris disks around planet-bearing stars is any different from that in debris disks 
without identified planets, nor that debris disks are more or less common (or more or less dusty) around 
stars harboring multiple planets compared to single-planet systems. Diverse dynamical histories may 
account for the lack of correlations. The data show a correlation between the presence of high-mass 
planets and stellar metallicity, but no correlation between the presence of low-mass planets or debris 
and stellar metallicity. Comparing the observed cumulative distribution of fractional luminosity to 
those expected from a Gaussian distribution in logarithmic scale, we find that a distribution centered 
on the Solar system’s value fits the data well, while one centered at 10 times this value can be rejected. 
This is of interest in the context of future terrestrial planet detection and characterization because it 
indicates that there are good prospects for finding a large number of debris disk systems (i.e. with 
evidence of harboring planetesimals, the building blocks of planets) with exozodiacal emission low 
enough to be appropriate targets for an ATLAST-ty^e mission to search for biosignatures. 

Keywords: infrared: stars — Solar system: interplanetary medium, Kuiper belt: general — stars: 
circumstellar matter, planetary systems, planet-disc interactions. 
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1. INTRODUGTION 

Planetesimals are the building blocks of planets, and 
mid- and far-infrared observations with Spitzer and Her¬ 
schel indicate that at least 10-25% of mature stars (ages 
of 10 Myr to 10 Gyr) with a wide range of masses (cor¬ 
responding to spectral types A-M) harbor planetesimal 
disks with disk sizes of tens to hundreds AU. This fre¬ 
quency is a lower limit because the surveys are limited by 
sensitivity. The evidence for planetesimals comes from 
the presence of infrared emission in excess of that ex¬ 
pected from the stellar photosphere, thought to arise 
from a circumstellar dust disk; because the lifetime of 
the dust grains (<1 Myr) is much shorter than the age 
of the star (>10 Myr), it is inferred that the dust can¬ 
not be primordial but must be the result of steady or 
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stochastic dust production generated by the collision, 
disruption, and/ or sub limatio n of p lanetesimals ( for re- 
views, see Wyat t 120081 Krivov [2010 Moro-Martin 2013 
Matthews et al. 2014p ' 


The Sun harbors such a debris disk, produced by the 
astero ids, comets, and Kuiper Belt objects (Jewitt et al. 


20091 with a dust production rate that has changed sig- 
nihcantly with time, being higher in the past when the 
asteroid and Kuiper belts ( Kbs) were more densely pop¬ 
ulated (Booth et al. 2009). Today, the Solar system’s 
debris disk is fainter than the faintest extrasolar debris 


disks we can o bserve with Herschel (Moro-Martin 2003 
Vitense et al. 2012), with a 3tT detection limit at 10- 


20 times the level of dust in the current KB (Eiroa et 


al. 2013 B. C. Matthews et al. 2015, in preparation). 
There is evidence of planetesimals around A- to M-type 
stars in both single- and multiple-star systems. These 
stars span several orders of magnitude difference in stel¬ 
lar luminosities, implying that planetesimal formation, a 
critical step in planet formation, is a robust process that 
can take place under a wide range of conditions. 

It is therefore not surprisin g that planets and debris 
disks coexist (Beichman et al. 

|2007b[ |2010[ M aldon ado et al. ' 


2005 

T 


Moro-Martin et al. 
Wyatt et al. |2012[ 
Marshall et al. 120141). However, based on Spitzer debris 


disk surveys, no statistical correlation has been found to 
date between the presence of known high-m ass planets 
an d debr is disks (Mor o-Ma rtin et al. 2007a Bryden et 
al. 120091 Kospal et al. 120091). These studies were focused 
on high-mass planets (>30 M 0 ) because, at the time, the 
population of low-mass planets was unknown. Overall, 
the lack of correlation was understood within the context 
that the conditions to form debris disks are more easily 
met than the conditions to form high-mass planets, in 
which case one would not expect a correlation based on 
formation conditions; this was also consistent with the 
studies that showed that there is a correlation between 
stellar metalli city a nd the presence o f mas sive planets 
(Santos et al. 2004 Fisher & Valenti 2005[ Maldonado 
et al. |2012 ), but there is no correlation between stellar 
me tallici ty and the pre sence of debris disks ( Greave s et 
al. 2006 Bryden et al. 2006 Maldonado et al. 2012). 

Recent results from the radial velocity surveys indicate 
that, similar to debris disks, there is no correlation be¬ 
tween the presence o f low- mass planets a nd ste llar metal¬ 
licity (Ghez zi et al. 2010 Mayor et al. [2011 Buchhave 
et al. 2012). This iriight indicate that the conditions to 
form low-rnass planets are more easily met than those 
to form high-mass planets. A natural question to ask is 
whether low-mass planets and debris disks are correlated. 

A correlation between terrestrial planets in the inner 
region of the planetary systems and cold debris dust has 
been predicted to exist based on a comprehensive set 
of dynamical simulations consisting of high-mass plan¬ 
ets, embryos, and inner and outer belts of planetesimals. 
These simulations find a strong correlation between the 
presence of cold dust and the occurrence of terrestrial 
planets because systems with cold dust imply a calm dy¬ 
namical evolution where the building blocks of low-mass 
planets have been able to grow and survive; on the other 
hand, systems with dynamically active high-mass planets 
tend to destroy both the outer dust-producing planetesi¬ 
mal belt and the building blo cks of the terrestrial planets 
(Raymond et al. 2011 2012). 


Herschel observations have opened a new parameter 
space that allows us to explore fainter and colder debris 
disks, improving our knowledge of debris disk frequency, 
in particular around later-type stars. In addition, since 
the Spitzer planet-debris disk correlation studies were 
carried out, a large number of low-mass planets have 
been detected, the frequency of which can now be char¬ 
acterized. Tentative detection of a correlation between 
low-mass p lanets and debris disks was presented in Wy¬ 
att et al. (2012) from a preliminary study based on a 
Herschel- DEBTUS subsample of the nearest 60 G-type 
stars, which was also seen in the volume-limited sample 
of radi al velo city planet host stars examined by Marshall 
et al. (2014). In this paper, we revisit the planet-debris 
disk correlation (or lack thereof) in the H ersche l DE¬ 
BRIS a nd D U NES s urveys (Matthews et al. 2010 Eiroa 
et al. (2010), 2013 B. C. Matthews et al. m prepara¬ 
tion) to assess whether the frequency and properties of 
debris disks around a control sample of stars are statis¬ 
tically different from those around stars with high-mass 
or low -mass planets. In a companion paper (Marshall 
et al. 2014), we describe the individual exoplanet host 
systenis, their debris disks, and the disk dependencies on 
planetary system properties such as planet semi-major 
axis and eccentricity. 

The selection criteria of the different samples used in 
this study are presented in Sectionj^ (with a discussion of 
biases in Section]^. A detailed discussion of the statis¬ 
tical analysis using Kolmogorov-Smirnov (K-S), Fisher’s 
exact, and survival analysis tests can be found in Section 

(regarding the frequency and properties of debris disks 
and their dependence on the presence of high-mass and 
low-mass planets), Section 0 (regarding the correlation 
with stellar metallicity), and Section [^ (regarding the 
distribution of the debris disk fractional luminosities). 
For a summary and discussion of our results the reader 
is directed to Section 0 

2. SAMPLE SELECTION 

Table 0 lists the selection criteria of the different sam¬ 
ples of stars used in our statistical analysis. Table [Ogives 
information on their stellar parameters, and Tablej^lists 
the observed fluxes and photospheric estimates at 100 
/rm and the strength of the excess emission. Detailed 
information on the procedures followed in this paper for 
source extraction, photosphere subtr action, and SEP fit¬ 
ting can be found in Kennedy et al. (2012a 2012b). 

All the stars included in this study are drawn from 
the Herschel DEBRIS and DUNES surveys. DEBRIS is 
an unbiased volume-limited survey for M-, K-, G-, F-, 
and A-type stars, where the volume limits are 8.6, 15.6. 


21.3, 23.6, and 45.5 pc, respectively (Phillips et al. 2010 
Oj B. G. Matthews et al. 2015, in preprara- 
iNEt 


Matthews 2010[ 
tion). The DUNES survey covers mid-F- to mid-K-type 
stars within 20 pc (irrespective of planet or debris disk 
presence), plus a handful of stars within 25 pc kno wn to 


harbo r planets and/or debris disks (Eiroa et al. 2010 
2 ^. 


2.1. Set 1: Control sample irrespective of planet and 
debris disk presence 

To maximize completeness from the DEBRIS and 
DUNES surveys, we selected for Set 1 all the FGK stars 
within 20 pc. 
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Table 1 

Sample description 


Set 


Description 


1 

2 

3 


4 

5 

6 
7 

ly, 2y, 3ay, 3by, 4y, 5y 

10, 2o, 3ao, 3bo, 4o, 5o 
toy, 2 oy 

loo, 2 oo 

11, 21, 3al, 3bl, 41, 51 
Ih, 21i, 3ah, 3bh, 4h, 5h 
It, 2t, 3at, 3bt, 4t, 5t 


FGK stars in DEBRIS and DUNES with distances <20 pc, ages > 100 Myr and no binary 
companions at <100 AU. 

Subset from Set 1 without known planets. 

Subset from Set 1 harboring high-mass planets with masses > 30 Mq. 

3a: for planets at > 0.1 AU. 

3b: for planets at < 0.1 AU. 

Subset from Set 1 harboring low-mass planets with masses < 30 M 0 . 

Subset from Set 1 harboring excess emission at 100 pm, i.e. with (.Fioo ~ -^*,ioo)/^iOO > 3). 
Subset from Set 1 with single planets. 

Subset from Setl with multiple planets. 

Subsets from Sets 1-5 with ages < 1 Gyr. 

Subsets from Sets l-5with ages > 1 Gyr. 

Subsets from Sets 1 and 2 with ages 0.1-5 Gyr. 

Subsets from Sets 1 and 2 with ages > 5 Gyr. 

Subsets from Sets 1-5 with metallicities smaller than the average [Fe/H] ^ -0.12. 

Subsets from Sets 1-5 with metallicities larger than the average [Fe/H] > -0.12. 

Subsets from Sets 1-5 with estimated dust temperature assuming a blackbody. 


The Spitzer surveys found that the upper envelope of 
the 70 fim debris disks emission shows a declin e over 
the ~ 100 Myr of a sta r’s lifetime (Bryd en et al. 2006 


Hillenbrand et al. 2008 Carpenter 2009). Therefore, to 
avoid introducing biases due to stellar age, we further 
restrict the control sample to stars with ages > 100 Myr 
(of the stars with known ages, only three were excluded 
because of y outh) Our stellar ages are o btained from 
Vican et al. (2012) and Eiroa et al. (2013). Stellar ages 
can be very uncertain, and individual systems may end 
up in the wrong age bin^^. However, for a statistical 
analysis such as the one in this paper, the best approach 
is to use an age database as ’’uniform” as possible. Our 
ages are based on gyrochronology, Ca// chromospheric 
emission (R’/fic) and X-ray flux, always in that order of 
priority, acknowledging the decreasing reliability of the 
corresponding age measu remen ts. Gyrochronology ages 
come from Vican et al. (2012) and are available for 17 
stars in our sample; they can be unreliable for young stars 
(< 300 Myr), but out of those 17 stars, only one star is 
in that age range. When several chromosphe ric age s are 
available, we favored the a ges in Eiroa et al. (2013) over 
those in Vican et al. (2012) because the latter were based 
on a literature search, whereas the former were derived 
using spectra obtained by the DUNES team and their 
innerly consistent estimates of Ca// activity index (out 
of the 162 c hromospheric ages used, 107 come from Eiroa 
et al. 2013). Stars without estimated ages were excluded 
from our analysis. 

We do not include A-type stars in this study because 
the planet searches around these targets are preferen¬ 
tially done around evolved A-type stars (classes III, III- 
IV, and IV) with low er jitt er and narrower absorption 
lines (Johnson et al. 2011), whereas the A-type stars 
targeted by DEBRIS are mam-sequence (class V). There¬ 
fore, we do not have information on planet presence for 
most A-type stars in the DEBRIS survey. Regarding M- 
type stars, 89 were observed by DEBRIS, three harboring 


C omparing, for example, t he ste llar ages in Sierchio et al. 
| |2014| l to those in Vican et al. ( | 2012 | l, among the 48 stars that 
these two studies have in common, we find that differences in 
ages are less than 50%, except for five stars: HD126660/HIP70497 
(80%), HD23754/HIP17651(83%), HD189245/HIP98470 (733%), 
HD20630/HIP15457 (70%) and HD101501/HIP56997 (84%). The 
age estimations are therefore broadly consistent 


planets, one of which also harbors a debris disk (GJ 581 - 
Lestrade et al. 2012). We do not include M-type stars in 
this study because of low number statistics and because 
they might probe into a different regime of planetesimal 
and planet formation than the FGK-type stars. 

The DEBRIS and DUNES surveys include single and 
binary/multiple stars. Previous studies indicate that 
there are differences in both disk frequency and planet 
frequency between singles and binaries, and these could 
introduce a bias in our statistical analysis . Reg arding 
disk frequency, Rodriguez & Zuckerman (2012) found 
that, out of a sample of 112 main-sequence debris disk 
stars, 25%±4% were binaries, significantly lower than the 
expected 50% for field stars, with a lack of binary systems 
at separations of 1-100 AU; for the debris disk hosts in 
the DEBRIS sample, the multiplicity frequency is ^ 28% 
(D .R .Rodriguez et al. 2015, in prep aration). Re garding 
planet frequency, Eggenberger et al. (2007 20111 carried 
out a survey with VLT/NAGO to look tor stellar compan- 
ions around 130 nearby solar-type stars and found that 
the difference in binarity fraction between the nonplanet 
hosts and the planet hosts is 13.2%±5.1% for binary sep- 
aratio ns < 100 AU. In a more recent study, Wang et al. 
(2014) compared the stellar multiplicity of held stars to 
that of a sample of 138 bright Kepler multiplanet can¬ 
didate systems, hnding also that, for the planet hosts, 
the binary fraction is signihcantly lower than held stars 
for binary semimajor axes <20 AU. An additional obser¬ 
vation is that, even within the giant planet regime, bi¬ 
naries tighter than 100 AU show a different distribution 
of masses, suggesting a different for matio n mechanism 
and/or dynamical history (Duchene 2010). In view of 
all these studies, we have excluded from our samples 96 
binary systems with semi-major axis <100 AU to avoid 
introducing a bias in our analysis. In doing that, we are 
nat urally excluding all circumbinary disks (Kennedy et 
al. 2012b D. R. Rodriguez et al. 2015, in preparation), 
limiting our analysis to those that are circumstellar. This 
seems appropriate because one would expect that the de¬ 
gree to which the dust is affected by planets (if present) is 
different, whether the dust is circumbinary or circumstel¬ 
lar, and this could again bias any potential planet-disk 
correlation. 

Differences in infrared background levels could intro- 
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duce a bias to the debris disk detection; however, both 
the D UNES and DEBRIS surveys excluded targets that 
were predicted to be in regions with high contamination 
from galactic cirrus^^. In addition, all the targets in Set 1 
have been inspected to exclude, to the best of our knowl¬ 
edge, sources subject to confusion. 

The total number of stars in Set 1 (FGK stars within 
20 pc, ages > 100 Myr and no binary companions at 
<100 AU) is 204. All the other star samples discussed 
in the subsections below are extracted from Set 1, i.e. 
they fulfill the same criteria with respect to stellar type, 
distance, age, absence of close binary companions, and 
nearby confusion. 

Table lists the planetary systems found within Set 1. 
There are 22 stars harboring planets and an additional 
three with unconfirmed planetary systems, namely HD 
22049 (e Eri), HD 10700 (r Get) and HD 189567. 

Even though the targets are located at a range of dis¬ 
tances (see Figure [^, we do not expect this to introduce 
a significant bias to the planet-debris disk correlation 
study presented in this paper for the following reasons. 
Regarding planet detection, the Doppler studies do not 
depend on distance (although their sensitivity depends 
on V magnitude and spectral type, and this may account 
for the closer distances of stars hosting low-mass planets 
only). Regarding debris disk detection, (a) the DUNES 
observations are designed to always reach the stellar pho¬ 
tosphere at 100 ^m to a uniform signal-to-noise ratio > 5; 
(b) we assess the planet-debris correlation using survival 
analysis that takes into account the upper limits from the 
DEBRIS survey; and (c) we use a distance-independent 
variable, the dust excess flux ratio (E))™ ~ -^star)/^star> 
where F™ is the observed flux at 100 ^m and Fgj.™ is 
the expected photospheric value at that wavelength. 


2 . 2 . Set 2: No-planet sample 

Set 2 is the subset of stars from Set 1 without known 
planets, as of August 2014. The number of stars in this 
set is 182 (179 if including the three unconfirmed plane¬ 
tary systems). 


2.3. Set 3: High-mass planet sample 

Set 3 is the subset of stars from Set 1 known as of Au¬ 
gust 2014 to harbor one or more planets with masses > 30 
M 0 (> 0.094 Mj„p). We call this the high-mass planet 
sample. The planetary system properties are listed in 
Table 1^ The number of stars in this set is 16 (17 if in¬ 
cluding the three unconfirmed planetary systems). Note 
that some of these systems also harbor low-mass planets. 
We chose this limiting planet mass because for stars har¬ 
boring planets > 30 M 0 , there is a correlation between 
th e prese nce of planets an d stel lar metallicity (Santos et 
al. 2004 Fisher & Valenti 2005). On the other hand, for 
stars harboring planets < 30 M 0 , there is no correlation 
between the pr esence of planets an d stel lar metallicity 
(Ghezzi et al. 2010 Mayor et al. 2011). This might 


indicate differences in the planet formation mechanism 
that may affect the planet-debris disk correlation. We 



Distance (pc) 


Figure 1. Distribution of distances. Top: Stars without known 
planets (Set 2). Middle: The line-filled colored histograms corre¬ 
spond to the high-mass planet sample (Set 3; in red, with hatching 
from the top-left to the bottom right), low-mass planet sample 
(Set 4; in green, with vertical hatching) and debris disk sample 
(Set 5; in blue, with hatching from the top-right to the bottom 
left). Bottom: Cumulative fraction of distances (same color code 
as above). 

further divide this set into two subsets: 3a (for planets 
with a > 0.1 AU) and 3b (for planets with a < 0.1 AU). 

2.4. Set 4- Low-mass planet sample 

Set 4 is the subset of stars from Set 1 known as of Au¬ 
gust 2014 to harbor one or more planets with masses < 
30 M 0 and no higher-mass planets. We call this the low- 
mass planet sample. There are six stars in this set (eight 
if including the three unconfirmed planetary systems). 

2.5. Set 5: Debris disk sample 

Due to the wavelength coverage of the DUNES and 
DEBRIS surveys^^, this study is focused on the 100 ^m 
emission. Set 5 is the subset of 29 stars from Set 1 with 
debris disks detected by Herschel at 100 fim, i.e. stars 
for which the signal to noise ratio of the excess emission 

plOO_ plOO 

is SNRdust > 3, where SNRdust = / °io°o 2 _^Too 7 > 

*^obs T^star 

and are the observed flux at 100 p,m and the 
estimated photospheric flux, respectively, whereas 
and (Ts™r their 1-cr uncertainties. The 70 /rm Spitzer 


The unconfirmed planethost star a Cen B was observed as 
part of the DUNES and Hi-Gal programs, but it was excluded 
from this analysis because its high background level does not fulfill 
the DUNES and DEBRIS selection criteria, and our analysis is 
intended to be unbiased. 


DEBRIS and DUNES utilized the simultaneous 100 pm and 
160 pm imaging mode as the basis for their survey data, with both 
teams taking additional data toward selected sources using the 70 
pm and 160 pm imaging mode of PACS and 250 pm, 350 pm and 
500 pm imaging with SPIRE as appropriate. 


































































5 


observations do not identify any additional debris disks 
within Set 1. This indicates that the 100 /rm emission is a 
good tracer of the cold KB-like dust, and we will use it as 
our reference wavelength. The analysis presented in this 
paper is limited to cold KB-like debris disks (where cold 
refers to debris disks detected at 70-100 ^m); we are not 
including the warm debris disks identified by Spitzer at 
24 /rm and with no excess at 70 /rm (under this category 
there is only one planet-bearing star, HD 69830). 

Note that there are several targets harboring debris 
disks and/or planets that were observed with Spitzer hvA 
were not observed by the Herschel DEBRIS and DUNES 
surveys because of their high level of background emis¬ 
sion. 


2.6. Sets 6 and 1: Single-/Multiple-planet sample 

Set 6 is the subset of stars from Set 1 known as of 
August 2014 to harbor single-planet systems, while Set 
7 is the subset of stars with multiple known planets. 


2.7. Sets ly-5y and lo-5o: Young/Old samples 

If debris disks evolve with time, and the samples com¬ 
pared have different age distributions, this will introduce 
a bias in our analysis. We therefore divide the samples 
into stars younger than 1 Gyr (labeled as Sets ly-5y) 
and stars older than 1 Gyr (Sets lo-5o; our sample has 
no hot Jupiters in Set 3o), limiting the comparison to 
sets of similar ages (i.e., within the o or y groups). We 
find that the distribution of ages in the samples consid¬ 
ered (Figure show that planet-bearing stars (Sets 3 
and 4) tend to be older on average than the stars in the 
no-planet sample (Set 2); this is because Gyr-old stars 
have low magnetic activity, implying lower levels of radial 
velocity jitter that facilitate the Doppler studies. While 
this might result in planet-bearing stars having fewer de¬ 
bris detections if debris levels decrease with age. Figure 

shows little evidence for evolution in disk detect abil ity 
with time, and this is discussed further in section 3.1 


2.8. Sets Ih and 11: High/Low metallicity samples 

To explore the role of stellar metallicity, we divide Set 
1 into two subsamples, a high-metallicity sample (Set Ih) 
and a low-metallicity sample (Set 11), using the midpoint 
of the metallicity distribution of Set 1, [Fe/H] = -0.12, 
as the dividing value. 


3. DEBRIS DISK FREQUENCY AND DUST FLUX RATIO 


The observed debris disk frequencies are listed in Ta¬ 
bles Hi and[^ Due to the small sample size, the statis¬ 
tical uncertainties are calculated using a binomial distri¬ 


bution rather than the v/ZV Poisson uncertainty (see the 
appendix of Burgasser et al. |2003|) . Table shows that 
the control sample (Set 1) has a debris disk frequency 
of 0.14j/Q;g3, similar to that found by the Spitzer surveys 


at 70 /im (Trilli ng et al. 2008 Hillenbrand et al. 2008 
Garpenter et al. |2009). 4’his result is also in agreement 
with Gaspar et al. (|2013|) who found a Spitzer incidence 
rate of 17.5% within the DUNES sample. 


3.1. Dependence on stellar age 

If debris disks evolve with time, and the samples com¬ 
pared have different age distributions within the decay 




Age (Gyr) 


Figure 2. Distribution of stellar ages. Top: Stars without known 
planets (Set 2). Middle: The line-filled colored histograms corre¬ 
spond to the high-mass planet sample (Set 3; in red, with hatching 
from the top-left to the bottom right), low-mass planet sample (Set 
4; in green, with vertical hatching) and debris disk sample (Set 5; 
in blue, with hatching from the top-right to the bottom left). Bot¬ 
tom: Cumulative distribution of stellar ages (same color code as 
above). 

timescale, this will introduce a bias in the comparison 
of the debris disk frequencies and dust flux ratios. As 
mentioned above. Figure [^indicates that planet-bearing 
stars (Sets 3 and 4) tend to be older on average than the 
stars in the control samples because they are preferen¬ 
tially targeted by the Doppler studies. 

To test for disk evolution, we divide the samples into 
stars with ages 0.1-1 Gyr (labeled as Sets ly-5y) and 
stars older than 1 Gyr (Sets lo-5o). We then compare 
the disk frequencies and dust flux ratios in the young and 
old samples. Set 2y and 2o (lines 9 and 14 in Table i). 
We do this exercise in the no-planet sample to minimize 
the effect of planet presence, as the goal is to check for 
disk evolution alone. Comparing Set 2y (with a disk fre¬ 
quency of 7/46 = 0.15) and Set 2o (with a disk frequency 
of 16/126 = 0.13) and using a binomial distribution, we 
find that detecting seven or more disks in Set 2y, when 
the expected detection rate is 0.13 (taking Set 2o as ref¬ 
erence, i.e. the expected number of disk detections is 
0.13-46) is a 39% probability event (24% if including the 
unconfirmed planetary systems - Table i lines 1 and 
2). This probability is not low enough to maim that the 
higher incidence rate in the young sample compared to 
the old sample is significant. 

The latter, however, does not take into account the 
uncertainty in the expected rate of the reference sample 
(in this case. Set 2o). The Fisher exact test is more 
appropriate in this regard. To carry out this test, we 
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Table 5 

Debris disk frequency (at 100 fim) 




Excluding unconfirmed planets^ 

Including unconfirmed planets®' 


Set 

No. of excesses 

Excess freq.^ 

No. of excesses 

Excess freq.'^ 







(at 100 fim) 


(at 100 /im) 

1 

1 

29/204 

^•^^+0.03 

29/204 

^•-^^+0.03 

2 

2 

24/182 

Q IQ—0.02 

'^•-^'^+0.03 

22/179 

'^•-^^+0.03 

3 

3a,b 

3/16 

0 ig“0.06 

^•^^+0.13 

4/17 

^•^'^+0.13 

4 

4 

2/6 

'^•'^'^+0.21 

3/8 

Q oy-0.13 
'^•'^'+0.18 

5 

5 

29/29 

29/29 

6 

6 

3/12 

'^•^^+0.15 

4/14 

'^•^^+0.14 

7 

7 

2/10 

^•^^+0.17 

3/11 

^•^'+0.16 

8 

ly 

7/48 

0 1 '5“0.04 

U-J-3+0.07 

7/48 

^•-*^^+0.07 

9 

2y 

7/46 

0 1 '5“0.04 
'^•-‘-^+0.07 

7/46 

0 1 ES“0.04 

^•-*-^+0.07 

10 

3aby 

0/2 

0 

0/2 

0 

11 

4y 

0/0 


0/0 


12 

5y 

7/7 


7/7 


13 

lo 

21/146 

0 14-0.02 
'^•-^^+0.03 

21/146 

0 14-0.02 
'^•-^^+0.03 

14 

2o 

16/126 

Q 10—0.02 

^•^'^+0.03 

14/123 

0 11“0.02 
^•-^-^+0.04 

15 

3abo 

3/14 

'^•^-‘-+0.14 

4/15 

^•^'+0.14 

16 

4o 

2/6 

^•'3'J+0.21 

3/8 

Q oy-0.13 
^•'^'+0.18 

17 

5o 

21/21 

21/21 

18 

6o 

3/10 

'^•'^^+0.17 

4/12 

Q oo-O.lO 

19 

7o 

2/10 

^•^^+0.17 

3/11 

'-'•^'+0.16 


^ Unconfirmed planetary systems are HD 22049 (e Eri), HD 10700 (r Get) and 
HD 189567. 

^ The statistical uncertainties are calculated using a binomial distribution. 


Table 6 

Dependence with stellar metallicity 


Set 

No. of stars 

No. with 

No. with 

No. with 


in set 

high-mass planets®' 

low-mass planets®' 

debris disks 



(> 30 M®) 

(< 30 M®) 

(at 100 /j.m) 

11 ([Fe/H] -0.12) 

61 

1 

3 (5) 

9 

Ih ([Fe/H] > -0.12) 

75 

14 (15) 

3 

17 


^ Excluding unconfirmed planetary systems around HD 22049 (e Eri), HD 10700 (r Get) and 
HD 189567. The parenthesis shows the result when including these three planetary systems. 


classify the stars in the two samples in two categories 
regarding disk presence: stars with disks (SNRdust > 3) 
and without disks (SNRdust < 3). The null hypothesis 
in this case is that both sets (2y and 2o) are equally 
likely to harbor disks. The test gives a 60% probability 
to find the observed arrangement of the data if the null 
hypothesis were true (Table lines 3 and 4). Note that 
the Fisher exact test can only reject the null hypothesis, 
never to prove it true. The Fisher exact test in this case 
does not identify any evolution in disk frequency within 
the timescale considered. 

A variable that is commonly used to characterize the 
strength of the disk emission is the dust flux ratio, 
(Ubs° - Uta?)/Uta?: where is the observed flux at 
100 fiui and is the expected photospheric value 

at that wavelength. Table lists the observed dust 
flux ratio for all the stars in our study. The Scr up¬ 
per limits (preceded by ”<” symbol) are given for stars 
without significant detected emission and are calculated 
assuming the observed flux is -I- Sct),™, for stars 


with 0 < F^^^/aZ 
< 0 . 


< 3, and 3 (Tq™, for stars with 


Figure^ shows the cumulative distribution of the dust 
flux ratio, whereas Figure shows its dependency with 
stellar age. To assess quantitatively whether the data 
show a decay with time, we carry out survival analysis. 
This is favored over the Kolmogorov-Smirnov (K-S) test 
because the latter does not deal with upper limits, and a 
significant number of the targeted stars have F’lno/^loo < 
3 (see Tablej^and down-faci ng arr ows in Figure®. Using 
ASURV 1.2^avalley et al. 19921, which impleme nts th e 
survival analysis methods of I'eigelson & Nelson (1985), 
we carried out the univariate, nonparametric two-sample 
Gehan, logrank, and Peto-Prentice tests to compute the 
probability that Sets ly and lo have been drawn from the 
same parent distribution with respect to the dust flux ra¬ 
tio. The results are listed in Table[^ line 5. The logrank 
test is more sensitive to differences at low values of the 
variable (i.e., near the upper limits), whereas the Gehan 
test is more sensitive to differences at the high end (i.e.. 
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Table 7 

Debris disk frequency (at 100 fim) as a function of spectral type 


Total^ F-type^ G-type^ K-type^ 

No of exceBsea ExCeSS freq.*> No of excesses b No. of excesses ExCeSS freq.‘> No of excesses ExCeSS freq.*’ 

iNo. oi stars , ^ JNo. or stars ^ JNo. or stars , No. or stars , 

(at 100 fim) (at 100 /i^m) (at 100 fim) (at 100 fim) 


1 

1 

29/204 

^•-^^+0.03 

10/46 

0 22”*^''^^ 
^•^^+0.07 

11/61 

n 1Q—0-04 
^•-‘-^+0.06 

8/97 

0.08;°°2 

2 

2 

24/182 

Q 10-0.02 

^•-^'^+0.03 

9/42 

'^•^-^+0.08 

7/48 

Q 1 C —0.04 
^•-*-^+0.06 

8/92 

0-09;S.o“4 

3 

3a,b 

3/16 

0 iq-0-06 

^•-^^+0.13 

1/4 

^•^^+0.25 

2/9 

^•^^+0.18 

0/3 

0 

4 

4 

2/6 

'-'•'J'J+0.21 

0/0 


2/4 

Q c—0.2 
'^•^+0.2 

0/2 

0 

5 

5 

29/29 

10/10 


11/11 

8/8 



6 

lo 

21/146 

0 14-0.02 
^•-^^+0.03 

8/33 

0 

^•^^+0.09 

7/49 


6/64 

u.uy_^0.05 

7 

2o 

16/126 

Q 10-0.02 

^•-^'^+0.03 

7/30 

'J-^'^+0.09 

3/37 

o-o87o:S6 

6/59 

'^•-‘-'^+0.03 

8 

3abo 

3/14 

^•^-*-+0.14 

1/3 

n 33“0.14 

u.00^0.29 

2/8 

^•^^+0.19 

0/3 

0 

9 

4o 

2/6 

Q OO-0.13 

0/0 


2/4 

f) K—0.2 
^•^+0.2 

0/2 

0 

10 

5o 

21/21 

8/8 


7/7 

6/6 



^ Excluding unconfirmed planetary systems around HD 22049 (e Eri), HD 10700 (r Get) and HD 189567. 
^ The statistical uncertainties are calculated using a binomial distribution. 


at the detections; Feigelson & Nelson 1985). The Peto- 
Prentice test is preferred when the upper liniits dominate 
and the sizes of the samples to be compared differ. The 
probabilities are not low enough to claim definitively that 
the two sets have been drawn from different distributions 
in terms of the dust flux ratio. However, given that they 
are in the 3-11% range to assess the role of planet pres¬ 
ence, we will take the conservative approach of limiting 
the comparison of disk frequencies and dust flux ratios 
to stars with ages > 1 Gyr (i.e., within Set lo). 


3.2. Dependence on planet presence 
3.2.1. High-mass planets 

To assess the effect of high-mass planets on the pres¬ 
ence of debris disks, we compare the disk frequencies in 
Set 3o (3/14 = 0.21) and Set 2o (16/126 = 0.13), limit¬ 
ing, for the reasons explained above, the comparison to 
the stars older than IGyr. Using a binomial distribu¬ 
tion, we find that detecting three or more disks in Set 
3o when the expected detection rate is 0.13 (taking Set 
2o as reference, i.e. the expected number of disk detec¬ 
tions is 0.13-14) is a 27% event; the probability drops to 
9% if including the unconfirmed planetary systems (Ta¬ 
ble]^ lines 9 and 10). Based on these numbers, there is 
no evidence that debris disks are more common around 
stars harboring high-mass planets compared to the aver¬ 
age population, in agreement with previous s tudies based 
on Spitze r obser vations (Mor o-Mart m et al. 2007a| Bry- 
den et al. 2009 Kospal et al. 2009). 

Classifying the stars in both samples (Sets 2o and 3o) 
into stars with and without disks and using the Fisher 
exact test, we find that there is a 41% probability to find 
the observed arrangement of the data if the null hypothe¬ 
sis were true, where the null hypothesis in this case is that 
the stars with at least one giant planet (Set 3o) and the 
stars without known-planet planets (Set 2o) are equally 
likely to harbor disks. This probability is 11% if includ¬ 
ing the unconfirmed planetary systems (Table ^ lines 
11 and 12). The Fisher exact test, therefore, does not 
identify any correlation between debris disk frequency 
and high-mass planet presence. To test how different the 
disk frequencies would have to be for a correlation to be 




Figure 3. Gumulative frequency of the dust flux ratio at 100 
/im. Top: only for the stars with significant detected emission 
(i.e., Fioo/c’’iOO > 3 - this panel is biased to large excesses because 
for stars with faint photospheres, they can be included only if the 
have large dust flux ratios). Bottom: for all the stars assuming an 
optimistic case, where the adopted flux ratio for the targets without 
significant detected emission is its corresponding upper limit, and 
a pessimistic case, where the adopted flux ratio is 0. Black is for 
the stars with ages > IGyr (Set lo) and red is for stars with ages 
< IGyr (Set ly). 


identified by the Fisher exact test, we carry out the test 
using Set 2o and a hypothetical Set 3o, varying in the 
latter the number of stars with and without disks: we 
find that the disk frequency for Set 3o would have to be 
about 2.8 times higher than in Set 2o. The identification 
of smaller differences in disk frequencies by the Fisher 
exact test is limited by low-number statistics. 

Using survival analysis, we address whether the dust 
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0.1 1.0 10.0 
Age (Gyr) 


Figure 4. Top: Dust flux ratio at 100 /im as a function of stel¬ 
lar age. The circles correspond to detections (i.e., -Fioo /<^100 > 
3), while the down-facing arrows correspond to upper limits (i.e., 
Fioo/<^(-^iOo) < 3). Black is for the stars without known planets 
(Set 2), red is for the high-mass planet sample (Set 3), and green 
is for the low-mass planet sample (Set 4). Unconfirmed plane¬ 
tary systems appear in orange. The larger open blue circles indi¬ 
cate which of those stars harbor excess emission at 100 fim (Set 
5). Bottom: Same as above but for the fractional luminosity, 
assuming a blackbody emission from the excess. The circles cor¬ 
respond to dust detections (i.e., stars with SNR^ust > 3, where 

plOO _pl00 

SNRnnst = obs star \ while the down-facing arrows corre- 
./^1002 , ^100 2 ^ 

V obs ■^‘^star 

spond to upper limits (i.e., SNR^ust < 3). 

flux ratio, Fj°g^/Fstar, is affected by the presence of high- 
mass planets. Figures and show the distribution of 
the dust flux ratio. The results from survival analysis 
(Table lines 15 and 16) indicate that there is a high 
probabmty that the high-mass planet sample (Set 3o) 
and the no-planet sample (Set 2o) have been drawn from 
the same population in terms of the dust flux ratio at 100 
fim (and the result holds if we include the unconfirmed 
planetary systems). The data do not show evidence that 
the disks around high-mass planet-bearing stars harbor 
more dust than those without known planets but with 
similar stellar characteristics. 

3.2.2. Low-mass planets 

We now repeat the exercise above for the low-mass 
planet sample, comparing the disk frequency in Set 4o 
(2/6 = 0.33) to that in Set 2o (16/126 = 0.13). Us¬ 
ing a binomial distribution, we find that detecting two 
or more disks in Set 4o when the expected number of 
disk detections is 0.13-6 (taking Set 2o as reference) is 
a 18% probability event; the probability drops to 5% if 
including the unconfirmed planetary systems (Table 
lines 17 and 18). Based on these numbers there is no 
firm evidence that debris disks are more common around 
stars harboring low-mass planets compared to the aver¬ 
age population. This test, however, does not take into 



Figure 5. Distribution of the excess flux ratio at 100 p.m for 
stars with significant detected emission (i.e., Fioo/cioo > 3). Top: 
The open black histogram corresponds to the stars without known 
planets (Set 2). Middle: The red filled histogram (with hatching 
from the top-left to the bottom right) corresponds to the high- 
mass planet sample (Set 3), while the green filled histogram (with 
vertical hatching) to the low-mass planet sample (Set 4). Bottom: 
Cumulative fraction (same color code as above). There are two 
stars outside the plotted range, one in Set 3a with = 

99.8 and another in Set 2 with F/OO/F/j.'lJ = 38.0. 


account the uncertainty in the expected rate of the ref¬ 
erence sample. 

The Fisher exact test gives in a 19% probability to 
find the observed arrangement of the data if the null 
hypothesis were true, where the null hypothesis is that 
the stars with low-mass planets only (Set 4o) and the 
stars without planets (Set 2o) are equally likely to harbor 
disks. The probability drops to 7% when including the 
unconfirmed planetary systems (Table lines 19 and 
17). We find that the disk frequency for Set 4o would 
have to be about four times higher than in Set 2o in 
order for the Fisher exact test to identify a correlation in 
our small subsample Set 4o. The identihcation of smaller 
differences in disk frequencies is limited by low-number 
statistics. 

The results from survival analysis (Table 1^- lines 23 
and 24) indicate that the probability that the low-mass 
planet sample (Set 4o) and the no-planet sample (Set 2o) 
have been drawn from the same population in terms of 
the dust flux ratio at 100 /rm is not low enough to claim a 
correlation (even when including the unconfirmed plan¬ 
etary systems). However, in this case, survival analy¬ 
sis might be unreliable because of the small sample size 
(N<10) of t he l ow-mass planet sample. 

In section[5]2]below we discuss that there are hints that 
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Figure 6. Cumulative frequency of the dust flux ratio at 100 
fim. Top: only for the stars with significant detected emission 
(i.e., -Fioo/crioo > 3). Bottom: for all the stars assuming an 
optimistic case, where the adopted flux ratio for the targets without 
significant detected emission is its corresponding upper limit, and 
a pessimistic case, where the adopted flux ratio is 0. Black is for 
the stars without known planets with ages > IGyr (Set 2o), red 
for stars harboring high-mass planets (Set 3o) and green for those 
harboring low-mass planets (Set 4o). The unconfirmed planetary 
systems are included under Set 2 (no-planet sample). 

the debris disk frequency around F-type stars might be 
higher than around G- and K-type, although this trend is 
not found to be statistically significant. However, given 
that none of t he F -type stars in our sample harbor plan¬ 
ets (see Figure [IT| because it is not possible to search to 
such low masses around them), to be conservative we now 
compare the low-mass planet sample to a control sam¬ 
ple that does not include F-type stars. We find that the 
binomial-derived probability that the disk frequencies of 
the low-mass planet sample and the no-planet sample 
(excluding the F’s) are similar is 9% (compared to 14% 
when including the F’s). The Fisher exact probability 
gives 12% (compared to 19% when including the F’s). 
Therefore, our conclusion that there is no evidence of 
correlation does not change when excluding F-type stars. 

In summary, our study does not show evidence of a 
correlation, but our conclusion is limited by the small 
sample size. 


with a disk frequency of 2/10=0.20) and using a bino¬ 
mial distribution, we find that detecting two or more 
disks in Set 7o when the expected detection rate is 0.30 
(taking Set 6o as reference, i.e. the expected number of 
disk detections is 0.30-10) is an 85% probability event 
(changing only slightly when including the unconfirmed 
planetary systems - Table lines 25 and 26). The data 
do not show any evidence that debris disks are more 
or less common around stars harboring multiple-planet 
systems compared to single-planet systems. The same 
conclusion results from the Fisher exact test (Table 
lines 27 and 28). Regarding the dust flux ratio, sur¬ 
vival survival analysis results (Table lines 29-34) indi¬ 
cate that the multiple-planet, single-jalanet and no-planet 
samples could have been drawn from the same popula¬ 
tion in terms of the dust flux ratio at 100 ^m (and the 
result holds if we include the unconfirmed planetary sys¬ 
tems). The data, again, do not show evidence of any 
correlation between planet multiplicity and the strength 
of the debris disk emission. 


3.2.4. Effect on the characteristic dust temperature 

We now assess whether there is any evidence that the 
debris disks around planet-bearing stars might be differ¬ 
ent from those around an average population of stars in 
terms of the characteristic dust temperature. Sets la¬ 
beled with a ”t” include only the stars with estimated 
dust temperatures (listed in Table [9). The calculation of 
the gra y-body dust temperatures is aescribed in Kennedy 
et al. (2012a) based on observations with a wide wave¬ 
length coverage. Figure shows the distribution of the 
characteristic dust temperature in the no-planet sample 
(Set 2t) and the planet samples (Sets 3t and 4t). The 
K-S test yields two values, D, a measure of the largest 
difference between the two cumulative distributions un¬ 
der consideration, and the probability of finding a D- 
value greater than the observed value; the latter is an 
estimate of the significance level of the observed value of 
I? as a disproof of the null hypothesis that the distribu¬ 
tions come from the same parent population; that is, a 
small probability implies that the distributions could be 
significantly different. The result from the K-S test is 
shown in Table (lines 35 and 36), showing a very high 
probability. The calculation of the probability is good if 
NiN 2 /(Ni-|-N 2 ) > 4, where Ni and N 2 are the number 
of stars in each set. However, if one wants to be conser¬ 
vative, it might be compromised when N < 20, as it is 
the case here. Based on the limited information we have 
so far, there is no evidence that the characteristic tem¬ 
perature of the debris disks around planet-bearing stars 
differs from the rest. 

4. CORRELATIONS WITH STELLAR METALLICITY 

Figure shows the distribution of stellar metallicity. 
To assess^e correlation with metallicity, we create Sets 
lm-5m, constituted by stars in Set s 1-5 with known 
met allicit ies^^ from Maldonado et al. (20121 and Eiroa et 
al. (2013). These sets are further divided into stars with 
high metallicities (Sets lh-5h) and those with low metal- 
licities (Sets 11-51), using the midpoint of the metallicity 


3.2.3. Planetary system multiplicity 

Comparing Set 6o (single-planet sample, with a disk 
frequency of 0.3) and Set 7o (multiple-planet sample. 


Regarding possible sou rces o f biases due to stellar age and 
distance, Maldonado et al. argued that because the stars 

are at close distances from the bun (in our case within 20 pc), it is 
unlikely that they have suffered different enrichment histories. 
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Figure 7. Distribution of the estimated black-body dust tem¬ 
perature for the stars with debris disk detections at 100 /im (i.e., 
SNRdust > 3). The open black histogram corresponds to stars 
without known planets (Set 2t), while the line-filled colored his¬ 
togram corresponds to stars harboring high-mass planets (Set 3t; 
in red, with hatching from the top-left to the bottom right) and 
stars harboring low-mass planets (Set 4t; in green, with vertical 
hatching). The top panel excludes unconfirmed planetary systems 
e Eri and r Get, while the bottom panel includes both planetary 
systems. 


Figure show the cumulative frequencies of the dust 
flux ratio and the fractional luminosity of Sets Ih and 11, 
showing that there is a dearth of debris disks with high 
dust flux ratios and high fractional luminosities around 
low-metallicity stars. However, the probabilities listed 
in Table (line 49) indicate that this trend is not sta¬ 
tistically significant. We cannot rule out the hypothe¬ 
sis that the high-metallicity and low-metallicity samples 
have been drawn from the same distribution in terms of 
the dust flux ratio. We conclude that the Fisher exact 
test and survival analysis do not allow us to identify any 
correlation between high stellar metallicity and debris 
disks. 

4.2. Planet presence 

Comparing the planet and no-planet samples in terms 
of stellar metallicity with the Fisher exact test (Table 
- lines 39-42), we find that in the case of giant planets, 
there is a 0.2% probability to find the observed arrange¬ 
ment if the stars without giant planets (Set Im-Set 3m) 
and the stars with giant planets (Set 3m) were equally 
likely to have metallicities > -0.12, whereas for low-mass 
planets (Set Im-Set 4 vs. Set 4) this probability is almost 
100% (the result holds when including the unconfirmed 
planetary systems). From the K-S test, the probability 
that the no-planet sample and the high-mass planet sam¬ 
ple could have been drawn from the same distribution in 
terms of stellar metallicity is 0.2%, whereas the probabil¬ 
ity that the no-planet sample has been drawn from the 
same distribution as the low-mass planet sample and the 
debris disk sample is much larger (49%; Table 1^- lines 
43-46). 

5. POSSIBLE BIASES INTRODUCED BY THE SAMPLE 
SELECTION 


distribution, [Fe/H] = -0.12, as the dividing value. Table 
lists how many stars are in each subset. 

4.1. Debris disk presence 

We now compare the debris disk frequencies in Set Ih 
(17/75 = 0.23) and Set 11 (9/61 = 0.15). Using a bino¬ 
mial distribution, finding 17 or more disk detections in 
Set Ih, when the expected detection rate is 0.15 (tak¬ 
ing Set 11 as reference, i.e. the expected number of disk 
detections is 0.15-75) is a 4% probability event (Table 
- line 37), indicating that the disk frequencies in the 
mgh- and low-metallicity samples might differ. This re¬ 
sult, however, does not take into account the uncertainty 
in the expected rate of the reference sample (in this case 
Set 11). From the Fisher exact test, we find that there 
is a 28% probability to find the observed arrangement of 
the data if the null hypothesis were true, where the null 
hypothesis in this case is that the stars without disks 
(Set Im-Set 5m) and the stars with disks (Set 5m) are 
equality likely to have metallicities > -0.12 (Table [^- 
line 38). From the K-S test, the probability that the no¬ 
planet sample (Set 2m) and the debris disk sample (Set 
5m) could have been drawn from the same distribution in 
terms of stellar metallicity is 33% (39% when including 
unconfirmed planetary systems; Table |^- lines 47-48). 

Regarding the strength of the excess emission, we use 
survival analysis to check if the low-metallicity and high- 
metallicity samples could have been drawn from the same 
population in terms of the dust flux ratio. Figures and 


5.1. Presence of undetected planets 

We now describe the potential biases that the sample 
selection could introduce in the statistical analysis de¬ 
scribed above. First, we assess whether the presence of 
unidentified planetary systems could affect our results. If 
we were to have many stars with high-mass planets in the 
control sample, Set 2, one could argue that a high-mass 
planet-debris disk correlation could have been present 
but hidden by all the “planet contaminants”. However, 
because the high-mass planet frequency is small, this 
seems unlikely. Due to the highe r freque ncy of low-mass 
planets (Mayor et al. 2009 2011 1 Batalha 2014 and refer¬ 
ences therein; Marcy et al. 20141 and references therein). 


we probably have many stars with low-mass planets in 
the control sample which have not been identified. This 
means that a low-mass planet-debris disk correlation may 
still be hidden in the data. We could avoid these biases 
by comparing the planet sets to a subset of stars in Set 
2 for which the presence of planets within a given period 
and mass has been ruled out by the radial velocity sur¬ 
veys. However, because nondetections are generally not 
made public by the planet search teams, the information 
to construct this no-planet stellar sample is not available. 


5.2. Distribution of spectral types 

By considering FGK stars to assess the planet-debris 
disk correlation, we are implicitly assuming that the disk 
frequency and the planet frequency do not differ signifi¬ 
cantly among these spectral types. 
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Figure 8. Distribution of stellar metallicities (logarithmic scale, 
with [Fe/H] = 0.0 for solar metallicity). The open black histograms 
correspond to stars without known planets and with known metal¬ 
licities (Set 2m). Top: The line-filled colored histograms corre¬ 
spond to stars harboring high-mass planets (Set 3m; in red, with 
hatching from the top-left to the bottom right), and low-mass plan¬ 
ets (Set 4m; in green, with vertical hatching). Middle: subset 
harboring excess emission at 100 pm (Set 5m; in blue, with hatch¬ 
ing from the bottom-left to the top-right). The stars with uncon¬ 
firmed planetary systems, e Eri and r Get, are included in Set 2m 
(no-planet sample). Bottom: Cumulative distributions of stellar 
metallicities (same color code as above). 


Tablej^and Figure [TT] show the distribution of spectral 
types in the samples under consideration. Let us limit 
the comparison to stars older than 1 Gyr (to avoid biases 
due to disk evolution), i.e. to the stars in Set lo (TableJ^ 
- line 6). For the F-stars, the disk fraction is 0.24 (8/o3 
disk detections): using a binomial distribution, finding 
eight or more disk detections, when the expected detec¬ 
tion rate is 0.14 (taking the disk frequency of the G-stars 
as reference, i.e. when the expected number of disk de¬ 
tections is 0.14-33) is an 8% probability event. While 
for the K-stars, with a disk fraction of 0.09 (6/64 disk 
detections), using a binomial distribution, finding six or 
more disk detections, when the expected detection rate 
is 0.14 is a 90% probability event. If we were to take the 
disk frequency of K-type as reference, for the F-stars, 
finding eight or more disk detections, when the expected 
detection rate is 0.09 (expected number of disk detec¬ 
tions of 0.09-33) would be a 0.8% event (Table [^- lines 
50-52). The latter seems to indicate there is a significant 
difference in disk frequencies between K-type and F-type 
stars. 

Eiroa et al. (2013) found that the frequency of disks in 




Figure 9. Cumulative frequency of the dust flux ratio at 100 
fim. Top: only for the stars with significant detected emission 
(i.e., Fioo/crioo > 3). Bottom: for all the stars assuming an 
optimistic case, where the adopted fiux ratio for the targets without 
significant detected emission is its corresponding upper limit, and 
a pessimistic case, where the adopted flux ratio is 0. Black is for 
the stars with metallicities larger than the average ([Fe/H] > -0.12; 
Set Ih) and red is for the stars with lower metallicities ([Fe/H] ^ 
-0.12; Set 11), independently of planet presence. 


the D UNES survey does not change significantly among 
FGK stars. The increased disk frequency for F-type 
stars found in our sample might have been biased to 
some degree by the shallower integration time of some 
of the DEBRIS targets, although the different Tes dis¬ 
tribution for the stars in the DEBRIS and DUNES sur¬ 
veys may also play a role (the former covering all FGK 
stars, whereas the latter covers mid-F to mid-K^®. Us¬ 
ing a larger sam ple oi Spitzer and Herschel observations, 
Sierchio et al. (20141 found no significant dependence 
with spectral type m the F4-K4 range. 

The test above does not consider the uncertainty in 
the expected rate of the reference sample. Glassifying 
the stars into those with and without debris disks and 
applying the Fisher exact test, we find that in this case 
the probability is not low enough to disprove the null 
hypothesis that the F-stars are equally likely to harbor 
disks as are the G-l-K stars (Table [^- lines 53 and 54). 

Regarding planet frequency, Doppler surveys indicate 
there is a correlation between high-mass planet frequency 
and spectral type that follows r oughly a linear increase 
with stellar mass (Johnson et al. 2010). From a compila- 


The spectral type dependence of the debris disk frequency 
within the DEBRIS sample will be studied in more detail by 
Sibthorpe et al. (in preparation). 
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Figure 10. Cumulative frequency of the dust fractional luminos¬ 
ity. Top: only for the stars with excess detections (i.e., stars with 
SNRdust > 3). Bottom: for all the stars assuming an optimistic 
case, where the adopted fractional luminosity for the targets with¬ 
out excess detections is its corresponding upper limit, and a pes¬ 
simistic case, where the adopted fractional luminosity is 0. Black 
is for the stars with metallicities larger than the average ([Fe/H] 
> -0.12; Set Ih) and red is for the stars with lower metallicities 
([Fe/H] ^ -0.12; Set 11), independently of planet presence. 


tion of Doppler surveys, Gaidos et al. (2013) suggest f(%) 
= -1.11 -I- 5.33 Mstar-/M0, for planets > S R® (masses 
> 95 M0 - see their Figure 8). For low-mass planets in 
the O.8-6R0 range, Kepler data indicate that among the 
FGK stars the planet frequency does not depend signif¬ 
icantly on the spectral type (Fressin et al. 20131. Ta¬ 
ble [7] and Figure 11 indicate that neither the high-mass 
nor low-mass planet frequencies within our sample reflect 
the above trends, with a higher incidence around G-type 
stars mostly likely because fewer F and Ks were searched 
for planets. This might skew slightly the disk incidence 
rate comparison for high-mass planets. Again, because 
nondetections are generally not made public, there is no 
way to circu mven t this issue. 

In Section [7(6| w e discuss how the conclusions change 
when excluding F-type stars from our analysis. 


6 . FRACTIONAL LUMINOSITIES AND COMPARISON TO 
THE SOLAR SYSTEM’S DEBRIS DISK 

Figure shows the cumulative frequency of the dust 
fractional luminosity. This variable is commonly used 
to characterize debris disk emission because it allows 
comparison of disks observed at different wavelengths; 
it is not very model-dependent as long as the wave¬ 
length coverage is good (as is the case in our samples). 
For stars with dust excess detections (SNRjust > 3), 
the fracti onal luminosity is calculated following Kennedy 
et al. (2012a 2012b). For stars with dust excess 


ages 


120 ■ Set 1 
100 - 
80 - 
60 ’ 

40 - 
20 ’ 
oL 



Spectral Type 


Figure 11. Distribution of spectral types for the different sets. 


nondetections (SNRdust < 3), the 3cr upper limit to 
the fractional luminosity is calculated from = 




, e^star —1 ^ 


following equation (4) in Be- 


ichman et al. (2006), and assuming the observed flux 
is Pohs + for stars with > 0^ and 

do-Q™, for stars with < 0. In this expression, 

X = where v is the frequency corresponding to 100 
fim, Tjtar = Teff is the effective stellar photospheric tem- 
perat ure, a nd Tdust is assumed to be 50 K (as in Eiroa 
et al. |2013|). 

The fractional luminosity can help place the debris disk 
observations in this study in the context of the Sola r 
system’s debris disk. Following Bryden et al. (2006), 
we compare the observed cumulative distribution of trac- 
tional luminosity to those expected from Gaussian dis¬ 
tributions in logarithmic scale, with average values of 
lOx, 3x, lx, and 0.1 x that of the Solar system’s debris 
disk, assuming for the latter a fractional luminosity of 
10-6.5_ avoid biases due to disk evolution, we limit 
the comparison to stars older than 1 Gyr (Set lo). The 
observed and Gaussian-derived cumulative distributions 


are shown in Figure 13 The bottom panel shows that 











































































































13 




Figure 12. Cumulative frequency of the dust fractional luminos¬ 
ity. Top: only for the stars with excess detections (i.e., stars with 
SNRjugt > 3). Bottom; for all the stars assuming an optimistic 
case, where the adopted fractional luminosity for the targets with¬ 
out excess detections is its corresponding upper limit, and a pes¬ 
simistic case, where the adopted fractional luminosity is 0. Black is 
for the stars without known planets with ages > 1 Gyr (Set 2o), red 
for stars harboring high-mass planets (Set 3o) and green for those 
harboring low-mass planets (Set 4o). The unconfirmed planetary 
systems are included under Set 2 (no-planet sample). 


the blue line exceeds the most optimistic case at low frac¬ 
tional luminosities. This means that we can reject the 
hypothesis that the median of the disk fractional lumi¬ 
nosity is 10 times that of the Sol ar syst em’s debris disk, 
in agreement with Bryden et al. (2006). The best fit to 
the data is a Gaussian centered on the Solar system value 
(magenta line in the top panel). This result is discussed 
in Section [7?7l 


7. SUMMARY AND DISCUSSION 

We have carried out a statistical study of an unbiased 
subsample of the Herschel DEBRIS and D UNES surveys, 
consisting of 204 FGK stars located at distances <20 pc, 
with ages > 100 Myr and with no binary companions 
at <100 AU. The main goal is to assess whether the 
frequency and properties of debris disks around a control 
sample of FGK stars are statistically different from those 
around stars with high-mass and low-mass planets. We 
find the following results. 


7.1. Disk evolution 

The Spitzer surveys found that the upper envelope of 
the 70 pjn debris disk emission shows a decline over the ~ 
100 Myr of a star’s lifetime, indicating that there might 
be a population of rapidly evolving disks that disperse 




Figure 13. Cumulative frequency of the fractional luminosity. 
The thick black histogram corresponds to the stars with ages > 
1 Gyr independently of planet presence (Set lo). Because we are 
interested in the cumulative frequency of the stars for fractional 
luminosities greater than the minimum observed value, in calcu¬ 
lating the cumulative distribution we adopt the pessimistic case, 
where the fractional luminosity for the stars without excess is 0. 
The blue, green, magenta and red lines correspond to theoretical 
distributions that assume a Gaussian distribution of fractional lu¬ 
minosities in logarithmic scale, with average values of lOx, 3x, lx 
and O.lx that of the Solar system, respectively, and assuming for 
the Solar system a fractional luminosity of 10“®’^. We fixed the 
cumulative frequency of disks with L^ust/L* > 10“^ at 10% ac¬ 
cording to the observed result (in set lo), implying 1-sigma widths 
for the theoretical distributions of 0.4, 0.8, 1.18 and 2.0, for the 
blue, green magenta and red lines respectively. Top: showing only 
the detected range; there are only three targets with fractional lu¬ 
minosities below 8T0~®, compromising the fit to the data in that 
low range because of low-number statistics. Bottom: the dot¬ 
ted line that coincides with the solid line corresponds to the the 
pessimistic case, where the adopted fractional luminosities for the 
targets without excess detections are taken to be 0, while the sec¬ 
ond dotted line on the upper part of the panel corresponds to the 
optimistic case, for which the upper limits are adopted. 


by 100 Myr. Our sample does not show clear evidence of 
disk evolution on the gigayear timescale. This is in agree¬ 
ment with the lack of disk evolution observed at 70 /rm 
in Spitzer surveys for stars older than 1 Gyr^^ (Trilling 
et ai. |2008[ Hillenbrand et al. |2008[ Carpenter et al. 


20091. In a recent study, using both Spitzer and Her- 
sctiel observations, and u sing a sample 2.5 times larger 
than ours, Sierchio et al. (2014) found that for disks with 
fractional luminosities smaller than 10“® there is a sig¬ 
nificant decrease in the debris disk frequency between 3 
and 5 Gyr. To look for evidence of disk evolution in the 


Compared to the 70 /rm observations, the 100 /rm emission 
might also trace dust located further out, where the collision times 
are longer; if this second population of dust exists, one would ex¬ 
pect even less evolution at this longer wavelength. 
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5 Gyr timescale that could bias our results, we have di¬ 
vided the sample into stars with ages 0.1-5 Gyr (labeled 
as Sets loy and 2oy) and stars older than 5 Gyr (Sets 
loo and 2oo). We then compare the disk frequencies and 
dust flux ratios in both subsamples (lines 6 and 8 in Table 

. The overall resulting probabilities are not low enough 
to claim that the two sets have been drawn from differ¬ 
ent distributions in terms of the dust flux ratio, nor that 
their disk incidence rates differ significantly. The Fisher 
eFisher exactxact test (line 7 in Table 1^ also indicates 
that both sets are equally likely to harbor disks. We 
therefore do not find evidence in our restricted sample of 
disk evolution in the 5 Gyr timescale. 

7.2. High-mass planet presence 

Our sample do not show evidence that debris disks are 
more common around stars harboring high-mass plan¬ 
ets compared to the average population. This is in 
agreement with the studies based on Spitzer observations 
that found no correlation between fractional luminosi¬ 
ties, Ldust/Lstar, an d the pr esence of high- mass p lanets 


(Moro-Martm et al. 2007a 


Bry den et al. 2009). Fig¬ 


ure 8 in Maldonado et al. (20121 also shows this trend. 


where the stars with discs and planets seem to be well 
mixed with stars with only disks in terms of the frac¬ 
tional luminosity, but they did not carry out any statis¬ 
tical analysis. This issue will be revisited using a larger 
sample that combines Herschel DEBRIS, DUNES, and 
SHARPS observations (G. Bryden et al. 2015, in prepa¬ 
ration) . 

Overall, the lack of observed correlation between high- 
mass planets and debris disks was understood within the 
context of the core accretion model for planet formation, 
where the conditions to form debris disks are more easily 
met than the conditions to form high-mass planets. This 
is in agreement with the metallicity studies that indi¬ 
cate that there is a correlation between high-mass planets 
and stellar metallicity, but no correlation between debris 
disks and stellar metallicity. Additionally, the presence 
of debris disks around stars with a very wide range of 
proper ties, fr om M-type (Kennedy et al. 2007 Lestrade 
et al . |201 2|) to the progenitors of white dwarfs (Jura 
2003 2007|, implies that planetesimal formation is a ro¬ 
bust process that can take place under a wide range of 
conditions. Therefore, based on formation conditions, if 
planetesimals can be common in systems with and with¬ 
out high-mass planets, there is no reason to expect a 
correlation between high-m ass planets and debris disks 
(Moro-Martm et al. 2007aI. 

Another factor contributing to the lack of a well- 
defined correlation with planet presence might be that 
the dynamical histories likely vary from system to sys¬ 
tem, and other stochastic effects need also to be taken 
into account, e.g. those produced by dynamical instabil¬ 
ities of multiple-planet syste ms clearing t he outer plan¬ 
etesimal belt (Raymond et al. 2011 2012), the planetes¬ 
imal belt itself trigge ring p lanet migration and instabil¬ 
ities (Tsiganis et al. 120051 Levison et al. 2011), or the 
stripping of planetesimals from disks during stellar fly¬ 
bys in the first 100 Myr, when sys tems a re still in their 
dense birth cluster (Lestrade et al. 2011). 

Another aspect that needs to be taken into account 
is that the planets detected by radial velocity surveys 
and the dust observed at 100 fim occupy well-separated 


regions of space, limiting the influence of the observed 
closer-in planets on the dust production rate of the outer 
planetesimal belt; there are long-range gravitational per¬ 
turbations produced by secular perturbations fr om sin gle 
planets on eccentric orbits (Mustill & W yatt |2009D or 
multiplanet systems (Moro-Martm et al. 2U07b| 2010) 
that allow close-in planets to excite outer planetesimal 
belts, but the timescale of the former may be longer than 
the age of the system, and the latter is limited to certain 
planet configurations. 

7.3. Low-mass planet presence 

In a preliminary study, and using a different subsam¬ 
ple of the Herschel DEBRIS survey, Wyatt et al. (2012) 
identified a tentative correlation between debris and the 
presence of planets with masses < 95 M m. Using a dif¬ 
ferent subsample, Marshall et al. (2014) also found evi¬ 
dence that stars with planets < 30 M^are more likely to 
harbor debris disks than are stars with planets > 30 Mq 
( 6/11 vs. 5/26). There are aspects related to the dynam¬ 
ical evolution of planetary systems that could result in a 
higher frequency of debris disks around stars with low- 
mass planets compa red to those with high-mass planets. 
Wyatt et al. (2012) discussed two alternative scenarios: 
(1) if the planets formed in the outer region and mi¬ 
grated inward, low-mass planets would have been ineffi¬ 
cient at accreting or ejecting planetesimals, leaving them 
on dynamically stable orbits over longer timescales; high- 
mass planets would have been more efficient at ejecting 
planetesimals, leaving behind a depleted population of 
dust-producing parent bodies. (2) Alternatively, if the 
planets formed in situ, the timescale for the planet to 
eject the planetesimals is shorter in systems with high- 
mass planets than with low-mass planets. However, the 
true migration histories of the systems studied may be 
significantly more complicated than the story portrayed 
under the two scenarios described above. For example, 
in our own Solar system, it is now well established that 
the ice giants, Uranus and Neptune, migrated outward 
over a significant distance to reach their current loca¬ 
tions, sculpting the trans- Neptu nian population as they 
did so (Hahn & Malhotra 20051) . 

In this paper we have used tne cleanest possible sam¬ 
ple of the Herschel DEBRIS and DUNES surveys to as¬ 
sess if the data at hand can confirm the tentative detec¬ 
tion of a low-mass planet-debris correlation. Contrary 
to the preliminary analyses mentioned above, here we 
have discarded stars without known ages, with ages < 
1 Gyr and with binary companions <100 AU, allowing 
us to rule out possible correlations due to effects other 
than planet presence. We find that the data do not show 
clear evidence that debris disks are more common around 
stars harboring low-mass planets compared to the aver¬ 
age population. However, having a clean sample comes at 
a price because the smaller sample size limits the strength 
of the statistical result: a positive detection of a corre¬ 
lation could have been detected by the Fisher exact test 
only if the disk frequency around low-mass planet stars 
were to be about four times higher than the control sam¬ 
ple. 

The planet-debris disk correlation studies can shed 
light on the formation and evolution of planetary systems 
and may perhaps help “predict” the presence of planets 
around stars with certain disk characteristics. Far from 
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being a closed issue, this correlation (or lack of) needs to 
be revisited. In the near future, G. Bryden et al. (2015, 
in preparation) will address this question using a sample 
that combines Herschel DEBRIS, DUNES, and SKARPS 
surveys, overcoming to some degree our limitations due 
to the small sample size. However, there are another two 
aspects that need to be improved upon and, with the 
data at hand, cannot be addressed at the moment: our 
ability to detect fainter disks and to detect or rule out 
the presence of lower-mass planets to greater distances. 

Regarding the disk detections, our knowledge of cir- 
cumstellar debris is limited: we only have detections for 
the top 20% of the dust distribution, assuming all stars 
have a remnant circumstellar disk at some level; limits 
closer to the KB-level are only possible for nearby F-|- 
type stars, and we are incapable of seeing exact analogues 
to our own Solar system leaving a large parameter space 
with no constraint on planet or dust properties. Future 
missions under consideration such as SPICA would im¬ 
prove things significantly: if its telescope is not descoped, 
the improvement in sensitivity would allow detection of 
photospheres not detected by Herschel, e.g. for M stars 
and for FGK stars at large distances; its noise would also 
be lower than Herschel, allowing it to detect fainter disks. 

Regarding the planet detection, the high frequency of 
low-mass planets indicates that we probably have many 
low-mass planet stars in the control sample which have 
not been identified, hindering our ability to detect a cor¬ 
relation. To overcome this problem, we rely on radial 
velocity surveys to gradually probe both to greater dis¬ 
tances and lower planet masses; but also critically impor¬ 
tant is that these teams make the nondetections publicly 
available so we can identify systems for which the pres¬ 
ence of planets of a given mass can be excluded out to a 
certain distance. 


not harboring planets, in which case it might not be ap¬ 
propriate to describe the dust excess emission with a sin¬ 
gle temperature. 


7.6. Stellar metallicity 

We find that there is no evidence that debris disks 
are more common around stars with high metallicities. 
Thi s is in agreement with previ ous studies (Greaves et 


al. 2006 Bryden et al. 2006). We find a dearth of 


debris disks with high dust mix ratios (also fractional lu¬ 
minosities) around low-metallicity stars, c onsist ent with 
the model of Wyatt, Glarke & Greaves (2007). How¬ 


ever, survival analysis tests indicate that this trend is not 
statistically significant and that we cannot rule out the 
hypothesis that the high-metallicity and low-metallicity 
samples have been drawn from the same distribution in 
terms of the dust flux ratio. 

The data confirm the well-known correlation between 
high metallicities and the presence of high-mass planets. 
On the contrary, we find no evidence of a correlation 
between high metallicities and the presence of low-mass 
planets. We therefore find the well-known positive corre¬ 
lation between the presence of planets and stellar metal¬ 
licity for stars with high-mass planets but no correlation 
for stars with low-mass planets only in agreem ent with 
extensive Doppler studies (Santos et al. 2004 Fi sher fc 

a Mayor et al. |MT1). 


Valenti 2005 Ghezzi et al 
Maldonado et al. 


( 2012 ) 


a larger stellar sam¬ 
ple and derived the metallicities in a uniform way. They 
found an increasing correlation with stellar metallicity 
from stars without planets and disks and stars with de¬ 
bris disks to stars with high-mass planets. They also 
concluded that the correlation with stellar metallicity is 
due to the presence of planets and not the presence of 
debris disks. 


7.4. Planetary system multiplicity 


Dynamical simulations by Raymond et al. (2011 2012) 
of multiple-planet systems with outer planetesimal belts 
indicate that there might be a correlation between the 
presence of multiple planets and debris. This is be¬ 
cause the presence of the former indicates a dynamically 
stable environment where dust-producing planetesimals 
may have survived for extended periods of time (as op¬ 
posed to single-planet systems that in the past may have 
experienced gravitational scattering events that resulted 
in the ejection of other planets and dust-producing plan¬ 
etesimals). It is of interest therefore to assess whether 
debris disks are correlated with planet multiplicity. 

Our sample does not show evidence that debris disks 
are more or less common, or more or less dusty, around 
stars harboring multiple-planet systems compared to 
single-planet systems. 


7.5. Dust temperature 

Based on the limited statistics, there is no evidence 
that the characteristic dust temperature of the debris 
disks around planet-bearing stars is any different from 
that in debris disks without identified planets. This is 
of course subject to detailed individual modeling, as the 
spatial dust disk distribution of the planet-bearing sys¬ 
tems might show more structural features due to gravita¬ 
tional perturbations compared to the disks around stars 


7.7. Fractional luminosity and comparison to the Solar 
system debris disk 

Comparing the observed cumulative distribution of 
fractional luminosity to those expected from a Gaussian 
distribution in logarithmic scale, we find that a distribu¬ 
tion centered on the Solar system value (taken as 10“® ®) 
fits the data well, whereas one centered at 10 times the 
Solar system’s debris disks can be rejected. 

This is of interest in the context of future prospects for 
terrestrial planet detection. Even though the Herschel 
observations presented in this study trace cold dust lo¬ 
cated at tens of AU from the star, for systems with dust 
at the Solar system level, the dust dynamics is dominated 
by Poynting-Robertson drag. This force makes the dust 
in the outer system drift into the terrestrial-planet re¬ 
gion. This warm dust can impede the future detection 
of terrestrial planets due to the contaminant exozodiacal 
emission, with its median level, its uncertainty, and shape 
of its distribution being some of the parameters that may 
affect the aperture size required for a telescope such as 
ATLAST to be able to char acteriz e biosignatures (see, 
e.g., Stark et al. 2014 Brown 2015). Ruling out a distri¬ 
bution of fractional luminosities centered at 10 times the 
Solar system level implies that there are a large num¬ 
ber of debris disk systems with dust levels in the KB 
region low enough not to become a significant source of 
contaminant exozodiacal emission. Comets and aster¬ 
oids located closer to the star are other sources of dust 
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that can contribute to the exozodiacal emission (and for 
those, Herschel observations do not provide constraints), 
but planetary systems with low KB dust-type of emis¬ 
sion likely imply low-populated outer belts leading to low 
cometary activity. These results, therefore, indicate that 
there are good prospects for finding a large number of de¬ 
bris disk systems (i.e. systems with evidence of harbor¬ 
ing planetesimals) with exozodiacal emission low enough 
to be appropriate targets for terrestrial planet searches. 
Dedicated warm dust surveys w ith th e Keck Interferome¬ 
ter Nuller ( Millan- Gabet et al. 2011), CHARA/F LUOR 
(Absil et al. [M^, VLTI/PIORTER (Ertel et al. [M4| 
and LBTl (under the HOSTS program) are shedding or 
will soon shed light on this issue. 

Even though the planetesimals detected by Herschel in 
the far infrared are located far from the terrestrial-planet 
region, their presence is favorable to the growth and sur¬ 
vival of terrestrial planets because these planetesimals 
indicate that the system has experienced a calm dynam¬ 
ical evolution, as opposed to an environment of dynam¬ 
ically active, high-mass planets. Such an environment 
would tend to destroy both the outer, dust-producing 
planetesimal belt and the planetesimals that might oth¬ 
erwise build the terrestrial pla nets. This c onclusion was 
the result of Raymond et al.’s (2011 2012) extensive dy¬ 
namical simulations consisting ot high-mass planet em¬ 
bryos and inner and outer belts of planetesimals. These 
simulations find that there is a strong correlation be¬ 
tween the presence of cold dust in the outer planetary 
system and the presence of terrestrial planets in the inner 
region, so a system with low levels of KB dust emission 
might also imply a dynamical history not amicable to 
terrestrial planets. The Solar system, in this case, would 
be an outlier, with a low-level of KB dust but a high 
number of terrestrial planets. It would be o f great inter¬ 
est to extend Raymond et al.’s (2011 2012) simulations 
to cover a wider range of initial conditions to further ex¬ 
plore this correlation, as it would enlighten the target 
selection for an ATLAST-type mission. 
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Table 2 

Stellar properties 


HIP 

HD 

GJ 

UNS 

Survey 

d 

(pc) 

SpT 

Teff 

(K) 

[Fe/H] 

(dex) 

Lbol 

(L©) 

Age 

(Gyr) 

Binary Sep.^ 

(AU) 

544 

166 

5 

G030A 

DUNES 

13.68 

G8 V 

5512 

0.15 

0.631 

0.17 


910 

693 

10 

F069A 

DUNES 

18.74 

F8 V Fe-0.8 CH-0.5 

6235 

-0.32 

3.280 

3.04 


1368 


14 

K115A 

DEBRIS 

14.99 

K7 

4089 

** 

0.109 

** 


1599 

1581 

17 

F005A 

DEBRIS 

8.59 

F9.5 V 

6013 

-0.29 

1.288 

4.03 


3093 

3651 

27 

K045A 

DUNES 

11.06 

KO V 

5261 

0.2 

0.515 

6.43 

A-B: 474 

3583 

4391 

1021 

G041A 

DEBRIS 

15.17 

G5 V Fe-0.8 

5845 

-0.2 

0.974 

0.84 

A-B: 252 



146.01 

G041B 

DEBRIS 

15.17 

G5 V Fe-0.8 

4228 

** 

0.019 

0.84 

A-B: 252 

3765 

4628 

33 

K016A 

DEBRIS 

7.45 

K2.5 V 

5006 

-0.24 

0.293 

5.31 


3821 

4614 

34A 

G004A 

DUNES 

5.95 

F9V - GOV 

5932 

-0.23 

1.209 

6.43 


3821 

4614 

34B 

G004B 

DUNES 

5.95 

F9V - GOV 

5932 

-0.23 

1.209 

6.43 


3909 

4813 

37 

F038A 

DUNES 

15.75 

F7 V 

6260 

-0.16 

1.770 

2.93 


4022 

4967 

40A 

K127A 

DEBRIS 

15.61 

K6 V k 

4136 

** 

0.109 

1.34 

A-B: 262 



40B 

K127B 

DEBRIS 

15.61 

K5 

2541 

** 

0.001 

1.34 

A-B: 262 

4148 

5133 

42 

K089A 

DUNES 

14.16 

K2.5 V (k) 

4913 

-0.16 

0.292 

3.64 


5862 

7570 

55 

F032A 

DEBRIS 

15.12 

F9 V Fe-l-0.4 

6144 

0.17 

1.950 

4.12 


7235 

9540 

59A 

G092A 

DEBRIS 

19.03 

G8.5 V 

5453 

0.01 

0.540 

2.20 


7513 

9826 

61 

F020A 

DUNES 

13.49 

F8V 

6224 

0.11 

3.446 

7.26 


7978 

10647 

3109 

F051A 

DUNES 

17.43 

F9 V 

6181 

-0.09 

1.571 

1.74 


8102 

10700 

71 

G002A 

GT 

3.65 

G8.5 V 

5421 

-0.43 

0.501 

5.82 


8768 

11507 

79 

K043A 

DUNES 

11.02 

K7 

3967 

** 

0.080 

0.67 


10798 

14412 

95 

G024A 

DUNES 

12.67 

G8 V 

5479 

-0.46 

0.431 

6.54 


11452 

15285 

98 


DUNES 

17.14 

K5V - M1.5V 

3921 


0.210 



11964 

16157 

103 


DUNES 

11.60 

K6Ve -MOVP 

3790 


0.092 

0.18 


12653 

17051 

108 

F046A 

DEBRIS 

17.17 

F9 V Fe-l-0.3 

6158 

0.07 

1.690 

1.52 


12777 

16895 

107A 

FOlOA 

DUNES 

11.13 

F7V - F8V 

6314 

0.03 

2.250 

7.92 


12777 

16895 

107B 

FOlOB 

DUNES 

11.13 

F7V - F8V 

6314 

0.03 

2.250 

7.92 


12843 

17206 

111 

F024A 

DEBRIS 

14.24 

F6 V 

6435 

0.04 

2.712 

0.71 


13375 


116 

K108A 

DEBRIS 

14.76 

K5 

4006 

** 

0.058 

** 


14445 

19305 

123 

K107A 

DEBRIS 

14.75 

K5 

3999 


0.101 

0.00 


15371 

20807 

138 

G018A 

DUNES 

12.02 

GO V 

5922 

-0.16 

1.009 

3.59 

A-B: 3717 

15457 

20630 

137 

GOllA 

DEBRIS 

9.14 

G5 V 

5738 

0.09 

0.846 

0.68 


15510 

20794 

139 

G005A 

DEBRIS 

6.04 

G8 V 

5500 

-0.34 

0.663 

6.56 


15799 

21175 

3222 


DUNES 

17.42 

KOV - KIV 

5087 

0.12 

0.535 

2.84 


15919 

21197 

141 

K122A 

DEBRIS 

15.39 

K4 V 

4534 

** 

0.220 

1.14 


16134 

21531 

142 

K061A 

DUNES 

12.51 

K6 V k 

4172 

-0.13 

0.126 

1.21 


16537 

22049 

144 

KOOIA 

GT 

3.22 

K2 V (k) 

5100 

-0.11 

0.337 

1.28 


16711 

22496 

146 

K079A 

DEBRIS 

13.59 

K5.0 

4194 

** 

0.121 

0.97 


16852 

22484 

147 

F022A 

DEBRIS 

13.98 

F8 V 

6031 

-0.07 

3.203 

7.59 


17420 

23356 


K087A 

DUNES 

13.95 

K2.5 V 

4982 

-0.12 

0.299 

7.35 


17439 

23484 

152 


DUNES 

16.03 

KIV -K2V 

5166 

0.05 

0.402 

0.76 




18 


Table 2 — Continued 


HIP 

HD 

GJ 

UNS 

Survey 

d 

(pc) 

SpT 

Teff 

(K) 

[Fe/H] 

(dex) 

Lbol 

(L©) 

Age 

(Gyr) 

Binary Sep.^ 

(AU) 

17651 

23754 

155 

F053A 

DEBRIS 

17.61 

F5 IV-V 

6646 

0.07 

5.158 

1.83 


18280 


156 

K124A 

DEBRIS 

15.52 

K7 

4121 

** 

0.114 

0.78 


19884 

27274 

167 

K067A 

DUNES 

13.05 

K4.5 V (k) 

4529 

0.06 

0.195 

0.00 


22263 

30495 

177 

G029A 

DUNES 

13.27 

G1.5 V CH-0.5 

5830 

0.04 

0.972 

0.55 


22449 

30652 

178 

F003A 

DEBRIS 

8.07 

F6V 

6538 

-0.01 

2.870 

3.04 


23311 

32147 

183 

K024A 

DUNES 

8.71 

K3+ V 

4755 

0.29 

0.283 

10.91 


23693 

33262 

189 

F012A 

DEBRIS 

11.64 

F6V 

6213 

-0.15 

1.496 

0.57 

A-B: 3743 

25544 

36435 

204.1 

G095A 

DEBRIS 

19.20 

G9 V 

5473 

0.06 

0.535 

0.40 


26394 

39091 

9189 

G085A 

DEBRIS 

18.32 

GO V 

6000 

0.1 

1.537 

6.04 


27072 

38393 

216A 

F006A 

DEBRIS 

8.93 

F6.5 V 

6374 

-0.09 

2.415 

3.45 

A-B: 860 


38392 

216B 

F006B 

DEBRIS 

8.93 

F9 V Fe-1.4 CH-0.7 

4905 


0.290 

3.45 

A-B: 860 

27188 


215 

K082A 

DEBRIS 

13.72 

K7 

4016 

** 

0.096 

** 


27887 

40307 


K065A 

DUNES 

13.00 

K2.5 V 

4844 

-0.45 

0.257 

2.60 


28103 

40136 

225 

F028 

DUNES 

14.88 

FOIV - F2V 

7000 

-0.11 

5.562 

0.64 


28442 

40887 

225.2 


DUNES 

16.39 

K5V - K6.5V 

4330 


0.290 



29568 

43162 

3389 

G056A 

DUNES 

16.72 

G6.5 V 

5582 

0 

0.711 

0.28 

A-B: 411 

32439 

46588 

240.1 

F056A 

DUNES 

17.88 

F8 V 

6244 

-0.03 

1.785 

7.05 


32480 

48682 

245 

F044A 

DUNES 

16.71 

F9 V 

6078 

0.16 

1.870 

6.32 


33277 

50692 

252 

G065A 

DUNES 

17.24 

GO V 

5924 

-0.11 

1.280 

5.45 


34017 

52711 

262 

G093A 

DUNES 

19.16 

GO V 

5915 

-0.14 

1.331 

5.54 


34052 

53680 

264 

G059C 

DUNES 

16.94 

GO V 

4314 


0.145 

6.24 

A-C: 3134 

34065 

53705 

264.1A 

G059A 

DUNES 

16.94 

GO V 

5826 

-0.39 

1.490 

6.24 

A-B: 357 

34069 

53706 

264.1B 

G059B 

DUNES 

16.94 

GO V 

5314 

-0.23 

0.492 

6.79 

A-B: 357 

35136 

55575 

1095 

F045A 

DUNES 

16.89 

F9 V 

5926 

-0.36 

1.497 

5.25 


37288 


281 

K099A 

DEBRIS 

14.53 

K7 

4054 

*= 4 ; 

0.073 

2.20 


37349 

61606 

282A 

K090A 

DEBRIS 

14.19 

K3- V 

4876 

0.06 

0.294 

0.46 

A-B: 827 



282B 

K090B 

DEBRIS 

14.19 

K3- V 

4074 

** 

0.102 

0.46 

A-B: 827 

38784 

62613 

290 

G062A 

DUNES 

17.18 

G8V 

5598 

-0.08 

0.629 

3.75 


40693 

69830 

302 

G022A 

DUNES 

12.49 

G8+ V 

5419 

0.14 

0.599 

6.43 


40702 

71243 

305 

F077A 

DEBRIS 

19.57 

F5 V Fe-0.8 

6702 

-0.25 

7.467 

0.25 


40843 

69897 

303 

F061A 

DUNES 

18.27 

F6 V 

6324 

-0.39 

2.507 

4.27 


42438 

72905 

311 

G036A 

DUNES 

14.36 

G1.5Vb 

5885 

-0.02 

0.999 

0.27 


42748 


319A&B 

K121AB 

DEBRIS 

15.37 

K7 

3778 

** 

0.083 

** 

A-B:1758 

43587 

75732 

324A 

K060A 

DUNES 

12.46 

KO IV-V 

5283 

0.19 

0.608 

8.43 

A-B: 1055 



324B 

K060B 

DUNES 

12.46 

KO IV-V 

2939 

** 

0.008 

8.43 

A-B: 1055 

43726 

76151 

327 

G068A 

DUNES 

17.39 

G3 V 

5815 

0.23 

1.028 

1.65 


44248 

76943 

332A&B 

F040AB 

DEBRIS 

16.06 

F5 IV-V 

6569 

0.08 

5.456 

0.83 

A-B: 455 

44722 


334 

K097A 

DEBRIS 

14.47 

K7 

3858 

** 

0.081 

1.00 


44897 

78366 

334.2 

G094A 

DUNES 

19.18 

GO IV-V 

6022 

0.03 

1.301 

1.20 


45343 

79210 

338A 

KOllA 

DEBRIS 

6.11 

K7 

3957 

** 

0.062 

0.70 

A-B: 106.7 

120005 

79211 

338B 

KOllB 

DEBRIS 

6.11 

K7 

3864 

-0.4 

0.068 

1.01 

A-B: 106.7 

46580 

82106 

349 

K064A 

DUNES 

12.89 

K3 V 

4781 

-0.44 

0.231 

0.61 


47592 

84117 

364 

F031A 

DUNES 

15.01 

F8 V 

6213 

-0.21 

2.018 

4.14 


48113 

84737 

368 

G086A 

DEBRIS 

18.34 

GO IV-V 

5968 

0.1 

2.629 

7.39 


49081 

86728 

376 

G040A 

DUNES 

15.05 

G4 V 

5796 

0.19 

1.369 

6.02 


49908 

88230 

380 

K005A 

DUNES 

4.87 

K5 

4090 

-0.16 

0.103 

6.60 


53721 

95128 

407 

G033A 

DUNES 

14.06 

GOV 

5921 

-0.03 

1.644 

4.93 


54646 

97101 

414A 

K056A 

DUNES 

11.96 

K7 V 

4468 

-1.5 

0.146 

4.11 

A-B: 408 



414B 

K056B 

DUNES 

11.96 

K7 V 

3483 

** 

0.051 

4.11 

A-B: 408 

55846 

99491 

429A 

G079A 

DEBRIS 

17.78 

G6/8 IH/IV 

5433 

*= 4 ; 

0.735 

3.79 

A-B: 504 

55848 

99492 

429B 

G079B 

DEBRIS 

17.78 

G6/8 IH/IV 

4942 

** 

0.325 

3.79 

A-B: 504 

56452 

100623 

432A 

K029A 

DUNES 

9.56 

KO- V 

5244 

-0.38 

0.364 

4.90 

A-B: 162.5 



432B 

K029B 

DUNES 

9.56 

KO- V 

9506 

** 

0.000 

4.90 

A-B: 162.5 

56997 

101501 

434 

G013A 

DEBRIS 

9.60 

G8 V 

5555 

0.01 

0.621 

1.66 


56998 

101581 

435 

K059A 

DEBRIS 

12.40 

K4.5 V (k) 

4557 

** 

0.152 

0.73 


57443 

102365 

442A 

G012A 

DEBRIS 

9.22 

G2 V 

5671 

-0.37 

0.839 

5.80 

A-B: 211 



442B 

G012B 

DEBRIS 

9.22 

G2 V 

2700 

** 

0.003 

5.80 

A-B: 211 

57507 

102438 

446 

G069A 

DUNES 

17.47 

G6 V 

5618 

-0.28 

0.699 

5.20 


57757 

102870 

449 

F009A 

DEBRIS 

10.93 

F9 V 

6130 

0.09 

3.594 

6.41 


57939 

103095 

451A 

K028A 

DUNES 

9.08 

K1 V Fe-1.5 

5167 

-1.12 

0.223 

4.61 


58345 

103932 

453 

K034A 

DUNES 

10.16 

K4-I- V 

4568 

0.16 

0.212 

0.44 


59199 

105452 

455.3 

F030A 

DEBRIS 

14.94 

FI V 

7053 

-0.13 

4.530 

0.68 


61094 


471 

K081A 

DEBRIS 

13.68 

K7 

3882 

** 

0.062 

** 


61174 

109085 

471.2 

F063A 

DEBRIS 

18.28 

F2 V 

6953 

-0.02 

5.092 

2.47 


61317 

109358 

475 

G007AB 

DEBRIS 

8.44 

GO V 

5929 

-0.2 

1.243 

3.96 


61901 

110315 

481 

K092A 

DEBRIS 

14.19 

K4.5 V 

4448 

** 

0.190 

6.66 


62207 

110897 

484 

F050A 

DUNES 

17.37 

F9 V Fe-0.3 

5939 

-0.53 

1.098 

6.24 


62523 

111395 

486.1 

G057A 

DUNES 

16.94 

G7 V 

5647 

0.22 

0.767 

0.37 


62687 

111631 

488 

K036A 

DEBRIS 

10.60 

K7 

4073 

** 

0.093 

0.68 


64394 

114710 

502 

GOlOA 

DEBRIS 

9.13 

GOV 

6072 

0.11 

1.418 

1.83 


64792 

115383 

504 

G073A 

DUNES 

17.56 

GOVs 

6043 

0.24 

2.187 

0.16 

A-B: 602 

64924 

115617 

506 

G008A 

DEBRIS 

8.55 

G7 V 

5597 

0 

0.836 

4.58 


65026 

115953 

508 


DUNES 

10.71 

KO, MOV - M3V 

3752 


0.120 

6.00 


65352 

116442 

3781A 

G050A 

DEBRIS 

15.79 

G9 V 

5248 

-0.14 

0.366 

7.81 

A-B: 416 

65355 

116443 

3782B 

G050B 

DEBRIS 

15.79 

G9 V 

5036 

-0.11 

0.289 

8.43 

A-B: 416 
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Table 2 — Continued 


HIP 

HD 

GJ 

UNS 

Survey 

d 

(pc) 

SpT 

Teff 

(K) 

[Fe/H] 

(dex) 

Lbol 

(Lo) 

Age 

(Gyr) 

Binary Sep.^ 

(AU) 

65721 

117176 

512.1 


DUNES 

17.99 

G4V - G5V 

5513 

-0.07 

2.989 

7.89 


66459 


519 

K041A 

DEBRIS 

10.93 

K5 

3986 

** 

0.063 

8.01 


66675 

118926 

521.1 

K109A 

DEBRIS 

14.76 

K5 

3876 

** 

0.066 

** 


67090 


525 

K070A 

DEBRIS 

13.19 

K5 

3852 

** 

0.047 

1.00 


67275 

120136 

527A 

F036A 

DUNES 

15.62 

F7V 

6826 

0.26 

3.018 

0.47 

A-B: 125.1 

67275 


527B 

F036B 

DUNES 

15.62 

F7V 

3580 

** 

0.032 

0.47 

A-B: 125.1 

67308 

120036 

1177 A 

K103A 

DEBRIS 

14.63 

K6.5 V (k) 

4116 

** 

0.115 

8.30 

A-B: 132.0 

67308 


1177B 

K103B 

DEBRIS 

14.63 

K6.5 V (k) 

4113 

** 

0.103 

8.30 

A-B: 132.0 

67487 

120467 

529 

K095A 

DEBRIS 

14.29 

K5.5 V k 

4293 

** 

0.169 

2.28 


67691 

234078 

532 

K091A 

DEBRIS 

14.19 

K5 

4169 


0.111 

6.36 


68184 

122064 


K032A 

DUNES 

10.06 

K3V 

4818 

0.1 

0.287 

3.34 


70218 


546 

K096A 

DEBRIS 

14.39 

K6 V 

4220 

** 

0.129 

0.40 


70319 

126053 

547 

G063A 

DUNES 

17.19 

G1.5 V 

5753 

-0.27 

0.820 

5.80 


70497 

126660 

549A 

F026A 

DEBRIS 

14.53 

F7V 

6328 

-0.08 

4.242 

1.08 

A-B:1005 



549B 

F026B 

DEBRIS 

14.53 

F7V 

3455 

** 

0.021 

1.08 

A-B: 1005 

71181 

128165 

556 

K072A 

DUNES 

13.22 

K3V 

4769 

-0.09 

0.240 

11.19 


71284 

128167 

557 

F039A 

DEBRIS 

15.83 

FSVwvar 

6889 

** 

3.557 

4.78 


71908 

128898 

560A 

AOlOA 

DUNES 

16.57 

A7p - FlVp 

7645 


11.263 




128898 

560B 

AOlOB 

DUNES 

16.57 







71957 

129502 

9491 

F062A 

DEBRIS 

18.28 

F2 V 

6759 

0.01 

7.474 

1.11 


73182 

131976 

570B 

K008B 

DUNES 

5.86 

K4 V 

3568 

-0.24 

0.069 

1.36 

A-BG: 141.8 

73184 

131977 

570A 

K008A 

DUNES 

5.86 

K4 V 

4607 

0.1 

0.210 

1.36 

A-BG: 141.8 

73996 

134083 

578 

F076A 

DUNES 

19.55 

F5 V 

6646 

0.05 

3.322 

3.95 


75277 

136923 


GIOIA 

DEBRIS 

19.60 

G9 V 

5369 

** 

0.497 

4.81 


76779 

139763 


K126A 

DEBRIS 

15.56 

K6 V k 

4161 


0.114 

2.93 


77257 

141004 

598 

G019A 

DEBRIS 

12.12 

GO IV-V 

5967 

0.06 

2.126 

5.49 


77760 

142373 

602 

G052A 

DEBRIS 

15.89 

GO V Fe-0.8... 

5897 

-0.5 

3.195 

7.39 


78072 

142860 

603 

FOllA 

DEBRIS 

11.26 

F7V 

6387 

-0.17 

3.097 

3.49 


78459 

143761 

606.2 

G064A 

DUNES 

17.24 

GO V 

5858 

-0.17 

1.780 

7.33 


78775 

144579 

611A 

K098A 

DUNES 

14.51 

KO V Fe-1.2 

5330 

-0.6 

0.439 

6.04 

A-B: 1020 



611B 

K098B 

DUNES 

14.51 

KO V Fe-1.2 

3372 

** 

0.003 

6.04 

A-B: 1020 

79248 

145675 

614 


DUNES 

17.57 

KOIV - KOV 

5336 

0.43 

0.653 

7.45 


79755 

147379 

617A 

K039A 

DEBRIS 

10.80 

K7 

4082 

** 

0.097 

1.01 

A-B: 697 

79762 


617B 

K039B 

DEBRIS 

10.80 

K7 

3345 

** 

0.029 

1.01 

A-B: 697 

80725 

148653 



DUNES 

19.66 

KIV - K2V 

5040 


0.636 

1.89 


82003 

151288 

638 

K031A 

DEBRIS 

9.81 

K5 

4418 

*= 4 ; 

0.125 

7.67 


83389 

154345 

651 

G088A 

DUNES 

18.58 

G8V 

5485 

-0.06 

0.619 

3.84 


83990 

154577 

656 

K080A 

DEBRIS 

13.63 

K2.5 V (k) 

4920 

-0.64 

0.221 

8.09 


84862 

157214 

672 

G035A 

DUNES 

14.33 

GOV 

5776 

-0.41 

1.281 

5.91 


85235 

158633 

675 

K062A 

DUNES 

12.80 

KOV 

5334 

-0.44 

0.417 

5.31 


85295 

157881 

673 

K021A 

DUNES 

7.70 

K5 

4201 

-0.03 

0.115 

2.20 


86796 

160691 

691 

G047 

DUNES 

15.51 

G3V/VI - G5V 

5787 

0.29 

1.821 

>1 


88601 

165341 

702A 

K007A 

DUNES 

5.08 

KOV 

5312 

0.05 

0.594 

1.05 


88601 

165341 

702A 

K007B 

DUNES 

5.08 

KOV 

5312 

0.05 

0.594 

1.05 


88972 

166620 

706 

K044A 

DUNES 

11.02 

K2 V 

5047 

-0.07 

0.346 

6.43 


89042 

165499 

705.1 

G075A 

DUNES 

17.62 

GO V 

5953 

0.09 

1.715 

5.40 


89211 

166348 

707 

K068A 

DEBRIS 

13.13 

K6 V (k) 

4225 

** 

0.126 

6.04 


91009 

234677 

719 


DUNES 

16.35 

K4V - K6V 

4200 

0.05 

0.232 

0.12 


92043 

173667 

725.2 

F073 

DUNES 

19.21 

F5V - F7IV 

6431 

0.04 

6.141 

4.74 


95149 

181321 

755 

G091A 

DEBRIS 

18.83 

G1 V 

5793 

-0.21 

0.791 

0.15 


95995 

184467 

762.1 


DUNES 

16.96 

KIV - 2V 

5027 

-0.22 

0.682 

7.48 


96100 

185144 

764 

G003 

DUNES 

5.75 

GOV - KOV 

5276 

-0.18 

0.427 

3.67 


97944 

188088 

770 


DUNES 

14.05 

K2V/VI - K4V 

4687 

-0.2 

0.770 

0.49 


98698 

190007 

775 

K063A 

DEBRIS 

12.86 

K4 V (k) 

4555 

** 

0.219 

0.59 


98959 

189567 

776 

G077A 

DUNES 

17.73 

G2 V 

5764 

-0.22 

1.024 

4.06 


99240 

190248 

780 


DUNES 

6.11 

G8IV 

5597 

0.3 

1.246 

8.30 


99701 

191849 

784 

K012A 

DEBRIS 

6.20 

K7.0 

3881 

** 

0.060 

0.85 


99825 

192310 

785 

K027A 

DEBRIS 

8.91 

K2-I- V 

5096 

0.09 

0.406 

7.50 


100925 

194640 

790 

G098A 

DEBRIS 

19.52 

G8 V 

5574 

0.06 

0.766 

4.78 


101955 

196795 

795 


DUNES 

16.72 

K5V - MO/IV 

4181 


0.331 

1.20 


101997 

196761 

796 

G037A 

DUNES 

14.38 

G8 V 

5486 

-0.3 

0.540 

5.63 


102186 

196877 

798 

K057A 

DEBRIS 

12.15 

K5.0 

4167 

** 

0.083 

0.49 


102485 

197692 

805 

F027A 

DEBRIS 

14.68 

F5 V 

6640 

0.03 

3.907 

0.31 


104092 

200779 

818 

K116A 

DEBRIS 

15.10 

K6 V 

4310 

** 

0.167 

1.12 


104214 

201091 

820A 

K002A 

DUNES 

3.49 

K5V 

4394 

-0.25 

0.144 

6.18 


104217 

201092 

820B 

K002B 

DUNES 

3.50 

K7V - MOV 

4002 

-0.39 

0.092 

8.45 


104440 


818.1C 

F079C 

DEBRIS 

19.81 

F9.5 V 

3370 

** 

0.008 

0.69 

AB-G: 142.6 

105090 

202560 

825 

K004A 

DEBRIS 

3.95 

K7.0 

3912 

** 

0.072 

2.56 


105858 

203608 

827 

F007A 

DEBRIS 

9.26 

F9.5 V 

6213 

-0.84 

1.522 

0.57 


106696 

205390 

833 

KIOIA 

DUNES 

14.62 

K1.5 V 

5013 

-0.2 

0.305 

2.01 


107350 

206860 

836.7 

G080A 

DUNES 

17.88 

GO V CH-0.5 

5992 

-0.2 

1.128 

0.32 

A-B: 772 

107649 

207129 

838 

G053A 

DUNES 

16.00 

GO V Fe-l-0.4 

5969 

-0.06 

1.282 

6.98 

A-B: 880 

108870 

209100 

845 

K003A 

DUNES 

3.62 

K4 V (k) 

4672 

-0.06 

0.210 

0.90 

A-BG: 1456 

109422 

210302 

849.1 

F064A 

DUNES 

18.28 

F6 V 

6463 

0.09 

2.883 

3.53 


110443 

211970 

1267 

K076A 

DEBRIS 

13.53 

K5.0 

4020 

** 

0.085 

3.78 
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Table 2 — Continued 


HIP 

HD 

GJ 

UNS 

Survey 

d 

(pc) 

SpT 

Teff 

(K) 

[Fe/H] 

(dex) 

Lbol 

(Lo) 

Age 

(Gyr) 

Binary Sep.^ 

(AU) 

111960 

214749 

868 

K077A 

DEBRIS 

13.55 

K4.5 V k 

4480 

** 

0.182 

3.96 


112447 

215648 

872A 

F043A 

DEBRIS 

16.30 

F7V 

6188 

-0.2 

4.780 

7.91 

A-B: 192.3 



872B 

F043B 

DEBRIS 

16.30 

F7V 

3370 

** 

0.024 

7.91 

A-B: 192.3 

112774 

216133 

875 

K088A 

DEBRIS 

14.13 

K7 

3854 

** 

0.056 

0.55 


113283 

216803 

879 

K019A 

DEBRIS 

7.61 

K4+ V k 

4578 

** 

0.187 

0.38 

A-B; 53498 

113421 

217107 


G102AB 

DEBRIS 

19.86 

G8 IV-V 

5645 

0.27 

1.173 

8.09 


113576 

217357 

884 

K022A 

DUNES 

8.21 

K5 

4079 

-1.5 

0.091 

4.90 


114361 

218511 

1279 

K114A 

DEBRIS 

14.99 

K5.5 V (k) 

4369 

** 

0.153 

1.26 


116215 

221503 

898 

K112A 

DEBRIS 

14.90 

K5 

4214 

** 

0.132 

8.07 

A-BG: 5047 

116745 

222237 

902 

K052A 

DUNES 

11.41 

K3+ V 

4743 

-0.24 

0.221 

6.60 


116763 

222335 

902.1 

G087A 

DEBRIS 

18.58 

G9.5 V 

5285 

-0.14 

0.439 

2.31 


116771 

222368 

904 

F021A 

DEBRIS 

13.72 

F7V 

6227 

-0.13 

3.512 

6.22 


120005 

79211 

338B 

KOllB 

DEBRIS 

6.11 

K7 

3864 

-0.4 

0.068 

1.01 

A-B: 106.7 


26976 

166B&C 

K006BC 

DUNES 

4.98 

K0.5 V 

3283 

** 

0.008 

4.30 

AB-C: 413 


^ Projected binary separation from Rodriguez et al. in preparation. 



Table 3 

Observed fluxes and excess properties 


21 


P 


9-3 


O ^ 
r 

O w 
>-1 0 


’-I 0 

b 

O m 
'H 0 


o 2-^ 
P°-S 


o 


-(M g<M - -CN c^(M<M(MiMC 

CO ^^^^^^^^ CO 

1 —I -<—I 


(MlOiM(M(M c3(M - - -(MCNlM -(MlM -C^l 

^ ^ ^ ^ ^ ^ ci .'' .'' <^i 


lO 

(M" 


' (M 




.o^oo<;o^co‘^^®'^■^^'^o^^~CN^~^^'^^iO'^'^'^*^^■^*^oococoo'^^^LO<^^'—lO^^'—ioo’^.;o^o^coi—i’^i>oo-^'^ 

^i-HlOi-H'^rHOOWCO^rHC^^^^^^O^^OQOC^rHrHOorHQ^^.^OOO(jJ^;Qt^l^(^QC^^QrH^^Ojv5oOQ^l^OO 


iO ;£) ^ 

lO o lO o o o 

'“’ ' O ' ' ■ 


^ ^ J:r ic CO 

O O o O o o o 

^l>Woc0Ol>00°5^*>®’^ 

•csi •i-HrH-!5^i-H(N^’^i-HrHC-J 

'^V'^^VVVVVvvVVV 


§ H ^ H H H 

W l> W o CO o 


jiO^o ;oloco lO coiio 

^'©OO'^^OOOOOOOOOOOO^OOO'^O^^^OO 
OOl lOOl I I I I I I I I I I lOi I lOiOOOl I 

WWcMCDWWc^^fMOilOO^COtM^i-IOOOW^OCqC^WioWWWt^oO 

2??5'^^<^^'^P‘©'^'-'OlO(^^<^^"^GOlO)t>OOpCO(MOO^>^^^0(;Cl^> 

0^00 • - 1 —ll>».• • -CO -ocoo 

. .,-(,-1 . .OOCS|IO) 1 -I 1 —I(M'^'-IIM(M^ 1 -I -IMI-IIO •'^ • • ■ ^ rsi 

cnoovv”°°VVVVVVVVVVVV^VVV' 


CO ^ t- CO 

o g o o 

l> W CO CO 
1-1 p 05 lO 


' V' 


CO CO 

lO o ^ o 

• 05 


VV VV^ V' 


CO lO CO ^ s 
O O O lO p 

HHHgH 

Tt GO ^ w CO 

O O lO o p 

I-H CO 00 • P 

VV VV^ V 


co^§io^oo^^;;^Ooi^PcooogcogioKoo^p!g2?it-S^^iocopg^^So^-Hg5a5t::^iOioJ^ 


V^' 


' ^ —v 

> Y Y o o o 


O O O OQ 

. V I oq V ' 


/ p l> O P o 


I IM l-l ^ ^ ^ 

; P P 9 P P d 


Tt^lOi-lOO'^C^a50<M^^COOoSc0^lO^COl005g'^2t>g^uOlOCOa5goO^^C^i^22o5t^(MlOiOlOi-l^-IiMi-lt>NlcO^T^ 
C0^i-HC0C0iO®5l>C0<M^-C0'^‘^^’^r-;*^rC-*^^r*^^-:O--^i-IC0a5^-l>*-. ^COCMP(M-CO<MOl>05-OOa5-COCO^ 


V 


O O O 1 




! CO CO *• . '^ 05 * 


’(M"Tt^05coS(M(M2Sfc’^05iM^2R 


;, 1.N (.N ‘ 1—I CN k1 

?(M(MPP05C^I>*^' 


]l>CO^P'^Pl0^lOo6ppo6o5CN' 


;a5co^OinPooo5^io‘ 

:oooooocoi-H*-.05a5>^io' 

Jc^COlOlOC^3‘^^'?*^! 


[OtMO^OO^ogS 

PcOCOiOi-i^(MPP 

Pi-hcoo6Pc^i>^P 



CO 

1—1 

CO 

CO 

p 

CM 

P 

p 


p 

t- 

p 

p 

lO 

p 

lO 

p 

uo 

p 

o 

p 

p 

lO 

p 

uo 

p 

CM 

p 

CM 

CM 

p 

X 

05 

p 

CO 

p 



05 

p 

CO 

CM 

X 

X 

X 

p 

CO 

p 

CM 

p 

lO 

X 


N- 

p 

CO 

p 

X 

p 

o 

p 

X 

p 

ti- 

p 

lO 

p 

05 

p 

p 

X 

p 

1—1 

p 

X 

iO X 
p p 

1—1 

X 

X 

N- 

p 

CO 

p 

o 

p 

X 


o 

CO 

CO 

CO 

CO 

CO 

CO 

CO 


CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

1—1 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

1—1 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 


CO 

CO 

CO CO 

CO 

CO 

CO 

CO 

CO 

CO 


-H 

-H 

-H 

-H 

-H 

p 

-H 

-H 

p 

-H 

-H 

p 

-H 

-H 

-H 

p 

-H 

-H 

-H 

-H 

p 

-H 

-H 

-H 

-H 

-H 

-H 

p 

-H 

p 

-H 

-H 

p 

-H 

-H 

-H 

-H 

-H 

p 

-H 

-H 

p 

-H 

-H 

-H 

-H 

-H 

-HP 

p 

p 

p 

p 

p 

p 

-ct^ 

lO 


lO 

<x> 

p 

CM 

00 

lO 

p 

CM 


CO 

CO 

1—1 

CM 

oq 

t^ 

p 

X 

o 

o 

05 

t- 


p 

05 

o 

p 

lO 

X 

CM 

lO 

lO 

p 

05 

CO 

p 

o 

o 

o 

X 

X 

CO 

X 

X 

X 

1—1 

1—1 

X 

p 

CM 


X 

05 

p 

05 

05 


CM 

CM 

CM 

CM 

p 

CM 

X 

CO 

p 

X 

1—1 

05 

CM X 
X 05 

05 

p 

CM 

X 

p 

CO 

lO 





CM 

CO 

CO 

CO 

GO 

CO 

CO 

I> 

00 

1—1 

CO 

CM 

CO 


CO 

CO 

05 

CM 

00 

1^ 

X 


CD 


lO 

05 

05 

CM 

CO 

00 


CM 

CM 

CO 

CM 

CO 

X 

CO 

-cr 

1^ 

05 

CO 

CM 



CM 

CM 


CO 

00 

X 

CO 

1—1 

1> 

CO 

X 


X 

lO 

l> 

CM 

i6 


lO 

CM 

lO 

lO 

lO 

CM 

05 

X 

CM 


00 

oq 

P 

o 

1—1 

CM 

P 

lO 

CM 

o 

CO 

lO 

05 

CO 

CM 

X 

P 

CM 

O 

05 

X 

CM 

p 

X 

X 

X 

p 

CO 

CO 

X 

X 

t- 

oq 

P 

P 

CM 

X 

CM 

lO 

CM 

N- 

X 

cq 

CO 

o 

o 

o 

lO 

cq 

X 

o 

cq 

1—1 

CO 

t- 

X 

oq 

05 

X 

CM 

1—1 


05 

cq 

X 

CM 

CM 

CO 

t- 

o 

cq 

1—1 

CM o 
CM t> 

p 

CM 

CO 

CM 

P 

P 

CM 

P 

X 

p 


P 

1—1 


p 

P 

P 

X 

1—1 

r-i 

cd 

cd 

o 

p 

P 


P 

P 

1—1 


P 

X 

cd 

P 


CM 

CM 



id 

P 

p 

cd 

id 

p 

P 

X 

P 

p 

id 

P 

P 

1—1 

X cd 

p 

p 

p 

p 

P 

P 

P 

P 

P 

p 

P 

P 

p 

p 

p 

P 

P 

p 

P 

P 

P 

P 

p 

p 

P 

p 

P 

P 

P 

P 

P 

P 

P 

P 

p 

P 

P 

P 

P 

P 

P 

p 

P 

p 

p 

P 

P 

P 

p 

p 

P 

P 

P 

PP 

lO 

lO 

05 

CM 

CM 

CO 

CO 

CM 

O 

CO 

CO 

p 

X 

p 

X 

X 

CO 

lO 

X 

X 

p 


lO 

05 

O 


O 

CO 

CO 

lO 

X 

CO 

CM 


CO 

t- 

CM 

X 

lO 

CO 


CO 

CM 

X 

X- 

X 

CM 

lO 

X 

lO 

i> 

CM 

-cr 

o 

CO 

CM 

X 

CM 

CM 

X 

CO 

X 

X 

CM 

CM 

P 

CM 

CO 

X 

o 

<05 

P 

CO 

CO 

CO 


X 

lO 

05 

O 

X Tf 

p p 

X 

CO 

05 

1—1 

p 

X 

CO 

;il 


i6 

GO 

X 

CO 

1—1 

X 

CO 

1 

X 

X 

X 

CM 

cd 

P 

X 

id 

CM 

cd 

05 

X 

cd 

P 

p 

id 

id 

CM 

p 

cd 

CM 

1 

P 

p 

X 

CM 

p 

lO 

id 

P 

X 

cd 

o 

cd 

X 

id 

cd 

05 

CM 

P 

P 

cd 

X 

CO 

<05 

P 

X 

X 

p 

p 


P 

00 

lo cd 


(M 


^ccc^cffl^cmccmc^cc^c^^ 

0051010^'—'’~lCO'^'^OOl>J>’05ir^C^O'-HCSCO'^ 
C0C0'~'O'^'^'^’~lOOc00C|C^G0C005(MiOO'^CN 

oo’^oPPooooo'^’^ooPoooPo 

OP^Pi^OO«OOP33^2itifeOPPO^O 


0051010^'—'’~lC0'^'^00l>J>’05ir^C^O'-HCSC0'^ COOO'^OOl>00.'^^ ^.•i....^w->.N«.' 

---- s-y^l ^ s^_.,—V (MCOOl>C^OO 

i-H O O O o o 


<c <; m <; <; < <; < <; <<;■<<;<;<; 

COOO'^OOl>OOi—llO 0^1—li—I05CS|N‘ 
"^1-|1-H(M00^1-I0 (MCOOl>C^OO 
OOOO'^’^OOO 1—lOOOoO 

fe[jHfefe;4^000 . ' ■ 


<;<;<<c<;<<;<;<;<;m<;<oo 

CO'^1>'05CO'^OC|^»OCOCOOC|LO^ 

i0CNC0<MoCSi-h05 00OO00C0o 

O'^OOoOoOOoOOOr- 


,0^l>ti-?3j5c0^3l>^^<Ml0^i-lO. __ 

^-'^^2^C0^^C0OO^^gc0^t-t-0505 


05 ic 00 § § ^ ^ ; 
n^ O O g g 1 


(MlOCOl>l>00C0O5^S!^f5^ 

ioiomcoi>t>oooo'^22^‘^i—' 

i-li-l^i-li-li-li-li-|OpTHrHcSl (M<M 

CS| 05 (^q CS| CSI 


OO lO 05 
lO CO CO 
^ CO 


OO'^'^COO 
(M ^ i-H 1-1 CO 
CO CO CO GO 05 
^ ^ ^ 


iO 1—( 1—I 05 CM 
CO CM CM O CM 
t> 00 00 05 O 
CO CO CO CO 


^^^,^t>Ot>CMlOl>i-llOlOCO l0t>O'^i0l>i-l05C0'^C0'^'^ ^lOCMl>CMlO^COCM l>COt> 

^PSS'^00^001010<35050 00c005t>05c0'^05 00i000l0 l>05iO'^COC0050505 OCOOO 

_::^r%^‘^l>^'^CMi-lO00 00CM COOOCOl>^i-llOO'^'^CO'^t> CM'^COi-lCM'Nt^OCOCO COi-HOO 

P^SSoO'-H'^lOCOt-COCOt- 05000^1-I1-HCMCMCMCOCOCO l>OOCMC0CO05 00 00 OOO 

u-Jl i-HCMCMCMCMCMCMCMCMCMCMCMCM CMCOCOCOCOCOCOCOCO 


GO CM lO CO 00 CM 00 f 


•’t>'^05COl>l>OOCO'^CO'^lOt>05i 


llOl>00'^'^COCMOOCM^COCOlOCOO 


^ 00 1-H 
I> l> OO 00 
CM CM CM CM 



Table 3 

Observed fluxes and excess properties (cont.) 


22 


3 8 
Q I 

"O 


■n 


o 

O “2 

o s 

i-HT! 

fin 

T! 

O 'S 




o 


(^"CN^CSC^CNCSC^CNC^CNCNC^ af(M(MCN c^"CSCN(^^CS^(^^C^(^^(^^C^(^^C^^''(^^(^^^''(^^C^^(^^ 




Q^O'^OOioQ9^^i— 

^ O 1-H 1-H (-H Q *0 C'^ ,_l (^ Q 1^ ^ CO CO ^ 


co^Sco^^^o^q^oq-^OM^Ol^^^^^q^a^ 


c^ o 


OOgOOOOOOOOOOgOOgOOOgOOOOOOOOooOOOOOOOOOOOOOOOOOgOOgOOO 


b- O W lO 05 'cO i-i CO 'iO 'cO 'lO O 
COCO'^CSCOCOCOO'^CSICOtMCO 

• • . 

COi-H •i-l<MiM(MCO(M'^lO(M^, 

vv^vvvvvvvvvv 


coa5WoGoa5Wco<Mcocoi>'^;3:^2ga5ioi>a5'^'^<Mcot>oooocoo50^ioi-ia5coW^aiW(Mo_ij 
'-ia5^CO^O>C^^(MC00005l>'^qqiOO(MOOCOOO<MO^l>(MO«000'^l>COCOCS|^COCS|^'^lO'^ 
rHi-H^-COi-HC^ •05<HrH^CsicO'^'^Tfo6o6o6'^'^C5"^rHi-Ho6rH(0-Jcsii-Hi-HCsii-HCO -lOrH -i-HC^^ 

VV'^VVV^VVVVVVvvVVVVVVVVVVVVVVVVVVV^VV'^VVv 


3 CO ^ wv 

co> o 05 o 


goS^;i, ... . . . 

00?^°6^00^^0^0000 


C005C^(^ 

^ CM Q 00 

Vi V' 




05 

, ^ c^i lo J 

• CO ^ 

' : 00 , ■ 




idddoooddd^^^dgo-dt^ddodd^ 


05 

w^o'^ooOoO 


O 05 

(M 05 ' 

-: t-; r 


v^ V' 


• ^ ^ CO 

' V VV 


CO 00 

00 CO'^^O’COoo'^i— 

'O dd°«cb°^d' 


i CO 00 ^ 1-1 ^ t 
1 00 05 (M • 05 • ( 


CO CO 00 lO S2 05 IIJ 

l00 05 00 C0t>C0<M00 05i-lC0oOlO'-l^lOo0^^‘30i . . . 

|^qqN.qqqqqqq^^^q^q^qo6®^^^®^^co4^--.N 

.y^OOiHOOOOO^y\/y^^OyO\/O^OO^y^OyC 






_ q q r _ 
lO lO O CO o t> ' 


I o CO GO ; 


.coo*- 

j E2 '^i ^ ' 


I O CO , 


•: 05 05 00 00t>^OCO .C0''cr''tl0)'Nrt>COOC0i 


iCO''ct^''crcO'i^^iO"^^iOiOO(M 


i O CO T 


! <M ^ 5 1> CO CO 

00 o ^ lO I> 05 r 


lO ^ o 


-i 




' o o o o o o o " 


' o o o o o ' 


; o o o o o 


o o o o " 


, 05 l> 
i CO lO 

; 1 -H o 


o o o o ^ 


+1 

+1 

d 

-H 

+1 

d 

-H 

+1 

+1 

d 

+1 

+1 

d 

+1 

+1 

+1 

+1 

d 

+1 

+1 

d 

-H 

+1 

d 

-H 

+1 

+1 

d 

+1 

+1 

d 

-H 

+1 

+1 

-H 

+1 

+1 

d 

-H 

+1 

d 

-H 

+1 

+1 

-H 

+1 

+1 

d 

-H 

+1 

d 

-H 

+1 

+1 

d 

0 

lO 

00 

10 

CO 

lO 

l> 

CO 

lO 

CM 

1—1 

q 

00 

CM 

lO 

0 

CO 

X 

CO 

05 

lO 

-t 

0 

CO 

CM 

0 

0 

05 

CO 

lO 

05 

CO 

X 

X 

lO 

05 

CO 

CO 

0 

X 

l> 

q 

CM 

X 

q 

b- 

b-* 

X 

q 

b- 

CM 

0 

q 

lO 

CM 

0 

X 

-t 

CM 

CM 

X 

X 

0 

0 

lO 

b' 

0 

lO 

X 

05 

CO 

CO 

X 

q 

lO 

-t 

X 

X 

q 

CO 

CM 

<35 

q 

-t 

q 

CO 

0 

X 

lO 

X 

lO 

CO 

05 

q 

lO 

06 



CO 

CM 

d> 

1—1 


CJ5 

CM 


CM 


05 

CM 


06 

CM 

05 

CO 

CO 

b- 

CM 

CM 

CO 

CO 

1—1 

CM 

16 

X 

X 


iro 

CM 

lO 

1—1 

-t 

X 

lO 

CO 

cj 

CD 

't 

rH 


0 

CM 

CO 


05 

b-^ 

06 

CM 

CM 

CO 

1—1 

X 


CO 

lO 


GO 

CM 

CO 

t- 

lO 

0 

CM 

CO 

<35 

CM 

lO 

CM 

0 

05 

lO 

CO 

CM 

05 

0 

q 

't- 

0 

CO 

l> 

X 

on 

05 


CO 

t- 

lO 

CO 

q 

05 

X 

CO 

CO 

X 

CO 

0 

q 

lO 

X 

-t 

lO 

lO 

05 

CO 

<35 

CO 

05 

CO 

0 

b- 

h- 

X 

CO 

b- 

iro 

q 

lO 

1—1 

CO 

0 

0 

X 

0 

lO 

X 

X 

b-; 

b- 

lO 

X 

X 

10 

05 

0 

X 

't 

<35 

-t 

CO 

X 

1—1 

<35 

<35 

b- 

X 

05 

b- 

-t 

1—1 

CO 

CM 

q 

CM 

X 

X 

q 

0 

X 

X 

X 

q 


rH 

CM 

1—1 


CM 

CM 



CM 

i-H 

CM 

CO 


CM 

CM 


rH 


CM 


CM 


CO 

1—1 


rH 

rH 


rH 


CM 


rH 

CM 



rH 

CM 

rH 

CM 

1—1 

rH 


r-i 

CO 

CM 

rH 

1—1 

CM 

CM 

X 

rH 

CM 

CM 

+1 

+1 

d 

-H 

+1 

d 

-H 

+1 

+1 

d 

+1 

+1 

d 

+1 

+1 

+1 

+1 

d 

+1 

+1 

d 

-H 

+1 

d 

-H 

+1 

+1 

d 

+1 

+1 

d 

-H 

+1 

+1 

-H 

+1 

+1 

d 

-H 

+1 

d 

-H 

+1 

+1 

-H 

+1 

+1 

d 

-H 

+1 

d 

-H 

+1 

+1 

d 

lO 

c:j5 

I> 

q 

CO 

CO 

0 

CO 

0 

CM 

't 

CM 

CM 

CO 

q 

h- 

-t 

CM 

X 

CO 

05 

0 

lO 

CM 

lO> 

CO 

05 

CM 

0 

05 

X 

CM 

lO 

CM 

q 

0 

0 

CO 

CM 

in 

X 

lO 

05 

b- 

X 

X 

lO 

CO 

05 

b- 

lO 

1—1 

1—1 

b~ 

b~ 

CO 

1—1 

CM 

CO 

CO 

X 

10 

05 

05 

1—1 

0 

CO 

X 

b- 

0 

1—1 

<35 

b~ 

b- 

X 

0 

't 

lO 

CO 

CO 

b- 

X 

b- 

CM 

b-* 

05 

X 

CM 

lO 

10 

in 

CO 

CO 

06 

CO 

CO 

CO 


0 

1 

06 



CO 

CO 

CD 

CO 

CO 

tV 

CM 


CO 

CM 


GO 

10 

lO 

1—1 

X 

X 

CO 

<35 

05 


b^ 

X 

b-^ 

CM 

CM 

CM 

1—1 

lO 

CO 

iro 

CM 

X 

CD 

't 

rH 

CM 

06 

1—1 


05 

lO 

CO 

b-^ 

tV 

CM 

CM 

-t 

iro 

g «6 

CM 

X 

<35 

lO 

-t 

05 

CM 

CM 


CO>CO'^l^'^'^0505ioC5000C|<Mt-i—iCO^OOOO'^l>'^i—li—l-^i—(COOiOC0COCO05 05 05 05 C005CM(M05 05 00''CfO'—'C0K^'^^O^>C0O 

ioiO'^coo5ioioiOTf05 05 05co<Mt^cocoq^2^og505^^^^q^qco^^t0^i!^^;i:^i2li?!iiJ^Q’^22^52^^2?'^r’iQ22lii? 


l0i0'^C005iOlbl0Ttiq05 05 C0 0Qt^C0C0^C0C0qO05qi—11—iqc0G0'^OC0l0l0b~l>6^0ClrHLbi—li—lqo^^^'C0C000^tr05lOlOC01—I 
00000000000000000^000^000000000000000000000000000200000 


■< m, 


CQ M 

Sq-i^-iococo^^2oo^^05oooi-i<^^cM<^co;2;^22 22'^®'^^‘^'^^^S2£^^^'^^£2£2'^'^ri^^^i>r'*>‘^*^2'^c- 


(M (M 


; CO CO 


CO 


CO 

CO 


CO CO CO 


^ ^ ^ 'It 


'It't 


-t 


(MOOOCi(M,-lOlOCOiO 
C00000051—lOOOOt> 
1—iiOCOCOl>COl>l>iO 
COCOOOOiMCOCOCOlO) 
■'t'lt'^LOlOiOLOlOiO 


CO O CO l> lO 
^ CO 05 O 
CO GO (N GO 05 
CM 05 i-H 05 
CO CO CO t> 


Jio •'t 
^ 05 
CO CO 


COOi—lCOt>l>OOOOOi—1 1 —ICM; 

COi-li-HO^COCMCOCMO 05 055 

COCMCM 1 —l^l>t>CMi-li-l 
G0O5O5CM'it''tCO00lOl> 05 05, 

I>l>l>00GO<X)000OO5O5 0505^ 


00 O lO CM CM 
CO t> 05 CO iO 
^ 00 O 05 'It 
CM CM CO CO iO 

00000 


-t 


lO GO lO l> lO 1 —I O 
00 lO 1 —( 05 05 CO 1 —I 
O CO CO 00 CO CO t> 
05 05 O O 1-1 1-1 ''t 
O O ^ 1-1 1-1 i-H 1 ^ 


X 

05 

0 

b- 

b~ 

CM 

lO 

<05 

CO 

X 

<05 


X 

CM 

X 

X 

X 

b- 

CO 

X 

CM 

b- 

X 

lO 

0 

CM 

X 


X 

1—1 

CO 

CO 

X 

CM 

b- 

X 

X 

b- 

b- 

05 

lO 

05 


-t 

b- 


b- 

X 

t- 


CO 

X 

X 

b- 


lO 

CO 

CO 

X 

X 

-t 

X 

<05 

0 

-t 

X 


X 

CM 

-t 

CM 

05 

't 


X 

05 


X 

0 

CM 

't 

-t 


lO 

<35 

05 


0 

iO 

X 

-t 

05 

05 

b- 


0 

0 

CM 

X 

05 

lO 



CM 

0 

0 

0 

0 


CM 

X 

b- 

CO 

b- 

X 


b- 

lO 

b- 

CM 

b- 

X 

X 


lO 

lO 


0 

05 

b- 

CO 

X 

X 


<35 

05 


lO 

b- 

05 

X 


0 

1—1 

X 

05 

CM 

lO 

CO 

X 

<35 

CM 

CM 

X 

-t 


-t 

-t 

iO 

b- 

b- 

X 

0 

0 

0 

CM 

CM 

X 

X 

't 

-t 


lO 


CO 

b- 

X 

05 

05 

X 

-t 

lO 

lO 

CO 

CO 

CO 

b- 

t- 

b- 

b- 

X 

05 


1—1 



CM 

CM 

CM 


CM 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 


't' 

-t 



-t 


-t 

-t 


't 


nb 


-t 


't- 

lO 

lO 

lO 

lO 

lO 

lO 

10 

lO 

lO 

lO 

lO 

lO 

lO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 



23 


a 

o 


ft 
CO w 


'V 

CD 

> 


o 


^ 8 

Q 'S 

"O 


T3 


O 

^ ■*^ 

C 

4^ 

O “2 

O 3 

i-HT! 

T! 

O^ 


7;°s. 


o 




CO CO 


rH rH rH ^ rH rH rH CO rH rH rH ^''rH rH ^''rH rH rH rH rH t-H rH rH t-H rH - - - ........ 


, CO ,-1 ^ CO r 


I CO , 


I CO , 




OO^^i—lOOOO’^COCOO'^'^'^O'^®^*^^ 
0^'ftOOI>’"^'-HOO’"^c"" 


lOi—lOt-^COb-'^O 

^oooc^ftoLO'ido'^oi^'^'^ooi-H 


giOco-^gggioioo^coooft(Ngcoi>coSqSS^"^?5t^§S§q^ 

<HPP'^ooi-H’^P’^’^^^’^o'^'^‘^cNoaio^'^o^o<oirHi-lo 


CO COCOiO COiOiO®-^iO»iOLOiO<lOLOl><lO'^ N- 

0^000*^00000000000000^^0 

hShhhShhhhhhhhhhhhhhSSh 

i-|W^i0)l>Wo<MO001>(MOOOi-Hi-Hl>C0CN|WP-lo 

co‘5cNi-ico2co'^^'^cO'-HOi>0'^coa5C^'^'®^oo 


I-H ^ O 1-H ' 

v^ V VV' 


'vvvvvvvvvvvvvv°°' 


V 


C000l>i0»i00l>-I>l>0'^l>i0-^N‘00»0)0 COOOOOOIOOOO 

oooooooooooooooooooogoooooooooo 

000'^lOC^t>0<MOlOO(MO'^COOOOi-HiOOWoooOOt>00<MOO'^ 

(MCOOa5(MO(Ma5'^(M'^Ot>COO'^0'^I>(M^t>0(MOC!OXOOi-lO 

i-HTfi-li-HrHC0i-H'^l0Oi-H-^Oi-HrHcdl0<0il^ii-H'l0rHi-lc^rHi-lTtt>i-ll> 

vvvvvvvvvvvvvvvvvvvv”vvvvvvvvvv 


S2^r^*^^2S*^OioJ2‘300ooooqco: 

OC^^ • -O • • ‘CO. 

^ •r^'-ifM-OOCO-OOC^Oi-Hi-Hr^'-'f 

^‘"TVV^VVV^VVVVViV 


) • ft »o • f 

,! (M o 7^, “f; o ' 
V V ' V V ' 






i O O CO ' 


_qq'- 
; o d o ' 


; o o 1 -H o ‘£^ )5 ' 

^ O O lO o ^ 'H o " 
^ rH d 9 d d ^ 


' O o o 


11—lOI>l>GOl> GO 1—I (NGO 1—O 

OCOOoO'^CiOoOOiOCOoOOOOOOft’COft’C^t—(COlOOTJ^OCO■^OcDOft’^CO^-cv^ql^^0^c^Oft'^>■COOO'—^OlOLOOOCOOl01—1(0) 
O(MO^--O_;^--C0---^--O^-O'-H-O'-H^-<MOOO--O<X)OO •^•OOOiOOOOOGOi-iOtMi-l- 

■C^ICO*!—ll>'^*OCO’—lC^0<l*C^-*O_'c*'^_;*'-C0'^''*''-H^*<MO —’. 


VV^VVV^VVVVV^V^ 


V V 


i-H o o o 


vv° 


(MCOftl>GO(MO^O;:5I’^OC^(MGo9co2!I^'^'^^®500i-Ht-'7*^o9(M<M^^l>'^COOCOGOOoS'^O^COO(MlOCoSo 
_ - ^.-|Cqi>iocOOOl>®i'7*.'^I>0'^'-HOOOCOGO^O°^OOOl>000'7''^iO"^OO^OOiO)inOOOO^GO 

d'pi-Hdo6i-Hco9d2Sl2dd9o6d9drH^dftI9d^dd'-Hcoddcoo6;i^dddo6do69t>P 


CO O ^ _ 

^'^COO'^^f^dcOiM 
o lO ft ft ft o 


!:G^OOO^OOOt-iMCOCOOl>i-HOOlO)Oo[^Ol>pI^< 

'^^qqq^qqq^qqrHqqqqqqqq^‘-.( 

‘^^ftftft^dddodddddddd^'^'^^^' 


; o o o o 


l> O I> lO T 

! o o o o ' 


ooooooooo" 


■H-H^-H-H^-H-H-H^-H-H^-H-H-H-H-H-H^-H-H^-H-H-H-H-H-Hq-H-H^-H-H-Hq-H-Hq-H-Hq-H-H-H^-H-Hq-H-H-H-H-H 


ft ^ GO O GO 

^ ^ q q ^ * 
^^dcoco^codi 


^ CO IM <M O CO lO i 


•9j9rHio^^9.qqqi>t-qi>qq^(Mo'd9^coooi>^q'^ 

]ddd'^d"!fddco^ 


-7 

CM 

q 

d 

o 

o 

o 

1—1 

o 

q 

m 

m 

CO 

q 

O 

q 

CO 

o 

O 

CO 

CTJ 

o 

00 

o 

q 

o 

q 

O 

q 

CO 

o 

o 

q 

00 

CO 

o 

00 

CO 

q 

CO 

q 

CO 

m 

q 

m 

q 

CO 

q 

1—1 

CM 

CO 

lO 

^3 


CM 

CM 

o 

CM 

CM 

00 

00 


O 

q 

CO 

q 

q 

00 

00 

o 

00 

-7 

00 

05 

CM 

q 

CO 

O 

00 

q 

o 

q 

lO 

i> 

CO 

CM 

CM 

00 

q 

05 

q 

00 

CO 

lO 

O 

q 



ft 

CM 

CM 

CM 

1—1 

CM 

ft 

CO 

d 

ft 

d 

CM 

CM 

ft 

1—1 

d 

CO 

CM 

1—1 

CM 

d 

CO 

CM 

CM 

CM 

ft 

CM 

CM 

CM 

CM 

CM 

1—1 

ft 

d 

d 

CO 

CO 

d 

1—1 

d 

CO 

1^ 

CO 


CM 

CO 

CO 

ft 

1—1 

CM 

d 

CM 

CM 

-H 

-H 

d 

-H 

-H 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

-H 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

m 

00 

q 

o 

q 

1—1 

CM 

q 

q 

CO 

CM 

CO 

q 

o 

q 

o 

CM 

o 

GO 

o 

q 

m 

lO 

CM 

-7 

q 

00 

00 

o 

CM 

o 

CM 

CO 

CM 

CO 

lO 

CM 

m 

CO 

m 

1—1 

1—1 

q 

CO 

00 

CO 

q 

CO 

CO 

CM 

q 

Oi 

q 

lO 

00 

CO 

q 

m 

o 

m 

05 

CM 

05 

q 

05 

q 

-7 

CM 

m 

CM 

00 

q 

CO 

o 

00 

m 

q 

1—1 

CM 

CO 

CO 

q 

m 

CO 

CO 

q 

CO 

m 

CO 

CO 

q 

o 

q 

00 

CO 

lO 

CO 

CO 

CM 

m 

d 

9 

m 


1—1 

d 

CM 

d 

d 

CO 

1 

d 

d 

00 

CM 

CO 

CO 

CO 

d 

d 

1—1 

m 

CO 

d 

CO 

d 

lO 

CO 

d 

CM 

CM 

d 

m 

CO 

d 

CM 

d 

CO 

CO 

d 

CO 

d 

CO 

d 

d 

CM 

d 

CO 


CM 

1—1 

d 

CO 

d 

d 

d 

CO 

05 

CM 

d 

1—1 

1—1 


d 

d 

CM 

d 


<< < ffl 

CO 00 o o 
l> O lO lO 

o o o o 

oo oo 




i-HojqooqqLOi— i(nocooo(MooO(MOOooo'z'^*^ 7:?'—•Tl'qoo 


1—ICJJOCOOCOC^LQi—(CN'^COOOCNO 
■^Ot^COCOOOOOCOOOcsiCNIb-CO'— 

Oi—iOr~if-^i—It—(OO0000000000000'~i'~*000000 


wt—lOOO'—I'^OO 


t—I CO o O l> O ' 
0000000000000'~I'~*000000 


< m 

at at 

CO CO 

o o 


o 


H <C CQ 

oooooooEl^^:^^^Mi:7ft^-^>(^^co 
t-t lo ^^ftft 


^ftoio^ftoo^oot^ 

iO> JiO JiO lO fti lO lO 


00<MCO*^22'^^S 

OOO^’^'^t— 

COOO OO 


<<:<:<:<:<: <: cQ << <: - 

i-lOOOlO<Mt—l^l>t>-t^-lO)O0 
COOOGOCOO(MoOO"^l>0 ft 

oooooopdooooo d 

<7 1—I 

OOt-HOlMlOCO^ftftO -t-o -lO 
C0lOiOl>l>t>O^^O!5Ot-li£Gio 
COOOCDOOO‘ftft’I>Ol>.t--.CNj>. 

l> t-* !> !>■ 


CO t> CO (M CO O 
GOt-HJO-^-^N- 
CO O O ^ ^ 

lO lO lO o o t> 


CO 

CM 

o 


CO 

CO 

o 


l> oo CO O 
CO l> O lO o 

o o o o 

O CM CO O 

CM CO CM CM CM 

t-H CM t-H t-H t-H 


CM"^OCMlO^OlOOiOlOOOGOl>i-l'^OOCr5t> 

(75CMCMlOlOCMl01>a>t>l>OOOOOOOt-li-iai 

l>OOCOCOt>'^OOCMCMCOCO'^Ot-lCMCO'7^ 

''crTfiOlOlOiCC0OI>I>l>t>I>l>l>00OOO 

cooooooooooooooooivt^t^ 


lOt^GOOOCMOl>COCOCO'^COOi—lO 
OOOOOt>l>OOCMOOl>OOl> 
1—It—iooooioa:iOOOt>ocoGOi>io 
OOOOGOOOCr5t—li-lT^COOt-ICMCMCO'^ 
CMCMCMCMCMCOCOCOCOCO'^'^'7*'^'^ 


1-I"^G0 l>CM'^'Ot>Ol>OCMOlO 
OOOOO lOOOOOa5t-t>lOON-lOl> 
1-icMO cr5t-ii-icr5CMt>CMi>0'^i> 
1—It—It—( t—lCOCOCOiOOl>l>GOOOOO 


mo COOOmO'^COt—1^1—It—lOOGOt-l>t—I 
t>l> moO'^I>i-lCOOOO'^'^CM<0'?tI>COCM 

coco mCMCOmCM^OOOmcOCOm'^CO'O'^OCO 

mi> G0i-i'^-^i>G0t>ommmm;0'^c0t-i 
-^lommLOiomm^mm^^ocob-oo 

ft ft ft 1—1 t—I t-H 1—1 t—I 1—I t—I 1—1 1—1 1—1 1—1 t—I CM 1—1 t—I 


GOmCMmcOCTiOCMlOmmi-lt—iCMCMt-HCJiCOCTi 
"^iOCOCMOOOC?i«OCOOOOON-'^t-iO'^'^ 
CMl>t>I>OCOOOOCMCMI>CO'0'OOCMOOt-l 
OOOOCMCOCO'^lOmmoOOOGOOOt—ICMIO 
I>l>t-GOOOOOGOOOOOOOGOOOOOGOOOGOOOO 



Table 3 

Observed fluxes and excess properties (cont.) 
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“ Observed PACS flux at 100 /rm with 1-cr uncertainty 
^ Estimated photospheric prediction at 100 fim with 1-cr uncer- 
tainty (<Ta\°°r)- 

^ Stars with significant detected emission have F^^g/cr^^® > 3. 

^ Observed flux ratio (Fobs/-^s\ar) dust excess flux ratio 

where = Fl^ - F^OJ. In both cases, the 3<7 

upper limits (proceeded by ”<” symbol) are given for stars with¬ 
out significant detected emission and are calculated assuming the 
observed flux is for stars with 0 < -F^bs/^obs ^ 

and 3crlg°, jriooy^ioo ^ g 

® Fractional luminosity of the dust excess emission. For stars 
with excess detections (SNRdust > 3), th e fractional lum inos- 
ity is calculated following Kennedy et ah (|2012a| |2012b||. For 
stars with excess nondetections (SNRdust < 3j, the '6a upper 
limit to the fractional luminosity is calculated from = 

100_pl0Q 

following equation (4) in Beich- 


man et ah 120061, and assuming the observed flux is Fq^s 
3(T^P'^, for stars with > 0, and 3(T^P'^, for stars with 

obs ’ obs ' obs ’ obs ’ 

^obs/^obs <0- fFis expression, x = where u is the fre¬ 
quency correponding to 100 ^m, Tstar = Feff is the effective stellar 
photospheri c tem perature and Tdust is assumed to be 50 K (as in 
Eiroa et ah 

^ Signal-to-noise ratio of the excess emission, given by SNRdust = 

plOO rlOO 


^ For label information, see Table ^ Systems that may be sub¬ 
ject to confusion are labeled as ’’(corn.?)”. The ’’set” classification 
of the systems with unconfirmed planetary systems -namely, HIP 
16537 (= e Eri), HIP 8102 (= r Get) and HIP 98959— are indicated 
in parenthesis. 

^ Note that this upper limit is based on the nondetection at 100 
/im; this star, however, has an excess emission at 8-35 pm with an 
inferr ed fractional luminosity of Ldust/Lstar = 2-10“'^ (Lisse et ah 
[2007) 1. 
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Table 4 

Planetary system properties' 


HIP 

HD 

GJ 

UNS 

Planet 

Name 

Mpisin(i) 

(Mjup) 

a 

(AU) 

e 

f^dust 

(AU) 

Set 

Ref.b 

3093 

3651 

27 

K045A 

b 

0.229 

0.29 

0.60 


3a,6 

(1) 

7513 

9826 

61 

F020A 

b 

0.669 

0.06 

0.01 


3a,7 

(2) 

7513 

9826 

61 

F020A 

c 

1.919 

0.83 

0.22 


3a,7 

(2) 

7513 

9826 

61 

F020A 

d 

4.116 

2.52 

0.27 


3a,7 

(2) 

7978 

10647 

3109 

F051A 

b 

0.925 

2.02 

0.16 

40.3 ± 5.9 

3a,6 

(3) 

12653 

17051 

108 

F046A 

b 

2.047 

0.92 

0.14 


3a,6 

(3) 

15510 

20794 

139 

G005A 

b 

0.008 

0.12 

0.00 

16.5 ± 7.5 

4,7 

(4) 

15510 

20794 

139 

G005A 

c 

0.007 

0.20 

0.00 

16.5 ± 7.5 

4,7 

(4) 

15510 

20794 

139 

G005A 

d 

0.015 

0.35 

0.00 

16.5 ± 7.5 

4,7 

(4) 

26394 

39091 

9189 

G085A 

b 

10.088 

3.35 

0.64 

51.3 ± 30.2 

3a,6 

(3) 

27887 

40307 


K065A 

c 

0.021 

0.08 

0.00 


4,7 

(5) 

27887 

40307 


K065A 

d 

0.028 

0.13 

0.00 


4,7 

(5) 

27887 

40307 


K065A 

b 

0.013 

0.05 

0.00 


4,7 

(5) 

40693 

69830 

302 

G022A 

c 

0.037 

0.19 

0.13 


4,7 

(6) 

40693 

69830 

302 

G022A 

d 

0.056 

0.63 

0.07 


4,7 

(6) 

40693 

69830 

302 

G022A 

b 

0.032 

0.08 

0.10 


4,7 

(6) 

43587 

75732 

324A 

K060A 

e 

0.026 

0.02 

0.00 


3a,7 

(7) 

43587 

75732 

324A 

K060A 

f 

0.173 

0.77 

0.32 


3a,7 

(7) 

43587 

75732 

324A 

K060A 

b 

0.801 

0.11 

0.00 


3a,7 

(7) 

43587 

75732 

324A 

K060A 

d 

3.545 

5.47 

0.02 


3a,7 

(7) 

43587 

75732 

324A 

K060A 

c 

0.165 

0.24 

0.07 


3a,7 

(7) 

53721 

95128 

407 

G033A 

c 

0.546 

3.57 

0.10 


3a,7 

(8) 

53721 

95128 

407 

G033A 

b 

2.546 

2.10 

0.03 


3a,7 

(8) 

55848 

99492 

429B 

G079B 

b 

0.106 

0.12 

0.25 


3a,6 

(3) 

57443 

102365 

442A 

G012A 

b 

0.051 

0.46 

0.34 


4,6 

(9) 

64924 

115617 

506 

G008A 

b 

0.016 

0.05 

0.12 

15.9 ± 1.5 

4,7 

(10) 

64924 

115617 

506 

G008A 

c 

0.033 

0.22 

0.14 

15.9 ± 1.5 

4.7 

(10) 

64924 

115617 

506 

G008A 

d 

0.072 

0.47 

0.35 

15.9 ± 1.5 

4,7 

(10) 

65721 

117176 

512.1 


b 

7.461 

0.48 

0.40 

14.0 

3a,6 

(3) 

67275 

120136 

527A 

F036A 

b 

4.130 

0.05 

0.02 


3b,6 

(11) 

78459 

143761 

606.2 

G064A 

b 

1.064 

0.23 

0.06 


3a,6 

(3) 

79248 

145675 

614 


b 

5.215 

2.93 

0.37 


3a,6 

(12) 

83389 

154345 

651 

G088A 

b 

0.957 

4.21 

0.04 


3a,6 

(13) 

86796 

160691 

691 

G047 

b 

1.746 

1.53 

0.13 


3a,7 

(14) 

86796 

160691 

691 

G047 

e 

0.543 

0.94 

0.07 


3a,7 

(14) 

86796 

160691 

691 

G047 

c 

1.889 

5.34 

0.10 


3a,7 

(14) 

86796 

160691 

691 

G047 

d 

0.035 

0.09 

0.17 


3a,7 

(14) 

99825 

192310 

785 

K027A 

c 

0.074 

1.18 

0.32 


4,7 

(4) 

99825 

192310 

785 

K027A 

b 

0.053 

0.32 

0.13 


4,7 

(4) 

113421 

217107 


G102AB 

c 

2.615 

5.33 

0.52 


3a,7 

(2) 

113421 

217107 


G102AB 

b 

1.401 

0.08 

0.13 


3a,7 

(2) 

Unconfirmed planetary systems^: 









16537 

22049 

144 

KOOIA 

b 

1.054 

3.38 

0.25 

36.0 

3a,6 

(15) 

8102 

10700 

71 

G002A 

b 

0.0063 

0.105 

0.16 

8.5 

4,7 

(16) 

8102 

10700 

71 

G002A 

c 

0.0097 

0.195 

0.03 

8.5 

4,7 

(16) 

8102 

10700 

71 

G002A 

d 

0.011 

0.374 

0.08 

8.5 

4,7 

(16) 

8102 

10700 

71 

G002A 

e 

0.013 

0.552 

0.05 

8.5 

4,7 

(16) 

8102 

10700 

71 

G002A 

f 

0.02 

1.35 

0.03 

8.5 

4,7 

(16) 

98959 

189567 

776 

G077A 

b 

0.0316 

0.11 

0.23 


3a,6 

(17) 


^ Planetary system properties from http: //ex oplane ts.org 
^ Orbit references are: (1) Wittenmyer et al. ) |2009|; (2 ) Wright et al. 
( |2009|l; (3 ) Butler et al. l |2006|l; (4) Pepe et ah l|20^y^ (5) Mayor et 
al. 1|20(5^; {^6) Lovis et al" (|2006|l; Endl et al. ||2012j|; (8) Gregory 
et al. pOlOy; (9) T inney et al. l |201l| |; (10) Vogt et al | |2010| ); (11) 
Brogi et al. I 2012^ ; (12) Wittenmyer et al. | |2007[ l; (13) Wright et al. 
l|2008|l; (14) Pepe et al. l|2007|l; (15) Hatzes et al. (|2000|l; (16) Tuomi 
et al. l |2013| l; (17) Mayor et al. j201l[ 

^ Uncontirrned planetary systems are HD 22049 (e Eri), HD 10700 (r 
Get) and HD 189567. 






























Table 8 

Results from the statistical tests 
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^ Including unconfirmed planetary systems? ”N” if those stars are 
included in the no-planet sample Set 2. ”Y” if they are considered 
planet-hosts (i.e., they are included in Sets 3 or 4 and 6 or 7). 

^ and are the total number of stars in each set 

(detections and nondetections). The number in parenthesis 
and are the number of stars in each respective set with 

upper limits (i.e the number of stars with nondetections for which 



logrank, and Peto-Prentice tests, indicating the probability that 


Sets A and B have been drawn from the same population in terms 


of the variable under consideration. 


^ K-S test probability. This is the probability that the cumulative 
distributions of the variable under consideration in Sets A and B 


differ by more than the observed value D, where D is a measure 
of the largest difference between the two cumulative distributions. 
A small probability implies that the distributions could be 
significantly different. 

® Fisher exact test two-tail probability. 

^ Using Poisson statistics, this is the cumulative probability of 
finding x or more disk detections (where x is the number of disk 
detections in Set B), when the expected rate is that of Set A. 

^ Using a binomial distribution, this is the probability of finding x 
or more disk detections (where x is the number of disk detections 
observed in Set B), when the expected rate is that of Set A. 


^ The Fisher exact test is calculated by dividing the samples into 
two groups: a high metallicity with [Fe/H] > -0.12 and a low 
metallicity with [Fe/H] ^ -0.12. The result is the probability that 
Sets A and B are equally likely to have the same distribution of 
high vs. low [Fe/H]. 


^ The Fisher exact test is calculated by dividing the samples into 
two groups: debris disks hosts, with a signal-to-noise ratio of the 


excess emission SNRdust 



> 3, and nondebris 


disks hosts, with SNRdust < 3. The result is the probability that 


Sets A and B are equally likely to harbor debris disks. 



HIP 

544 

1368 

4148 

5862 

7978 

8102 

15510 

16537 

16852 

17420 

17439 

22263 

23693 

26394 

28103 
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Table 9 

Debris disks properties (detected at 100/rm) 


HD 

GJ 

UNS 

'T'cold a 

■pcold a 

(At) 

166 

5 

G030A 

86.2 ± 2.0 

8.3 ± 0.4 


14 

K115A 

29.0 ± 3.2 

30.5 ± 6.8 

5133 

42 

K089A 

29.2 ± 2.6 

49.0 ± 8.7 

7570 

55 

F032A 

73.8 ± 23.8 

19.8 ± 12.8 

10647 

3109 

F051A 

49.1 ± 3.6 

40.3 ± 5.9 

10700 

71 

G002A 

80.0 ± 

8.5 ± 

20794 

139 

G005A 

61.8 ± 14.1 

16.5 ± 7.5 

22049 

144 

KOOIA 

35.0 ± 5.0 

36.0 ± 

22484 

147 

F022A 

98.0 ± 7.7 

14.4 ± 2.3 

23356 


K087A 

59.3 ± 83.3 

12.0 ± 33.8 

23484 

152 


41.0 ± 

29.0 ± 

30495 

177 

G029A 

70.6 ± 2.7 

15.3 ± 1.2 

33262 

189 

F012A 

115.0± 11.7 

7.2 ± 1.5 

39091 

9189 

G085A 

43.3 ± 12.7 

51.3 ± 30.2 

40136 

225 

F028 

149.0± 

8.4 ± 


32480 

48682 

245 

F044A 

51.9 ± 3.1 

39.3 ± 4.8 

42438 

72905 

311 

G036A 

87.2 ± 9.5 

10.2 ± 2.2 

43726 

76151 

327 

G068A 

87.0 ± 19.6 

10.4 ± 4.7 

61174 

109085 

471.2 

F063A 

37.4 ± 1.9 

124.6± 13.4 

62207 

110897 

484 

F050A 

53.7 ± 8.3 

28.2 ± 8.8 

64924 

115617 

506 

G008A 

66.8 ± 3.1 

15.9 ± 1.5 

65721 

117176 

512.1 


lOO.Oib 

14.0 ± 

71181 

128165 

556 

K072A 

42.5 ± 59.7 

21.0 ± 59.1 

71284 

128167 

557 

F039A 

126.8± 34.1 

9.1 ± 4.9 

85235 

158633 

675 

K062A 

62.0 ± 16.2 

13.0 ± 6.8 

107350 

206860 

836.7 

G080A 

86.6 ± 8.7 

11.0 ± 2.2 

107649 

207129 

838 

G053A 

44.1 ± 1.6 

45.2 ± 3.4 

114361 

218511 

1279 

K114A 

30.6 ± 3.3 

32.4 ± 7.1 

116771 

222368 

904 

F021A 

51.3 ± 29.1 

55.1 ± 62.5 


^ ^dust stars with 100 pm excesses, calculated fol¬ 

lowing Kennedy et al. ^2012a[ |2012b[ using the full spectral energy 
distribution. 








